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laser  and  a  pulsed  XeF  laser.  Mode  formation  was  deduced  from  observation  of 
the  far-field.  The  analytical  effort  consisted  of  a  fast  Fourier  transform 
(FFT)  wave  optics  computation  which  modeled  the  far-field  pattern  evolution. 

Results  agree  with  the  temporal  mode  formation  geometrical  theory  of  j 

Anan'ev  when  laser  saturation  is  achieved  at  the  predicted  mode  formation  time  1 
in  a  CO2  laser.  Rapid  saturation  at  one-third  this  time  results  in  significant! 
earlier  mode  formation  compared  to  the  geometrical  theory  for  a  XeF  laser. 
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I.  INTRODUCTION  AND  HISTORICAL  DEVELOPMENT 


A.  Background  of  Laser  Transverse  Mode  Formation 

The  evolution  of  an  ideal  laser  beam  from  the  onset  of  spontaneous 
emission  to  the  production  of  steady  state  transverse  resonator  modes  charac¬ 
teristic  of  a  diffraction  limited  output  depends  upon  resonator  geometry  and 
possibly  upon  saturable  gain.  This  phenomenon  is  a  particular  interest  in  the 
discussion  of  short  pulse  high  gain  lasers.  In  these  lasers  a  sizeable  amount 
of  energy  is  often  lost  during  formation  of  the  lowest  order  transverse  mode, 
leading  to  poor  far-field  focusing. 

Historically,  this  problem  devolves  from  two  kinds  of  laser  resonators, 
the  stable  and  the  unstable.  The  first  kind  of  resonator  to  be  used,  the 
conditionally  stable  resonator,  was  developed  as  an  extension  of  the  Fabry- 
Perot,  or  planar  mirror  resonator.  Later,  the  concave-planar  configuration, 
now  in  widespread  use,  was  developed.  Many  lasers  at  that  time  had  features 
which  conformed  to  the  following  parameters:  continuous  wave  (CW)  operation, 
relatively  small  beam  diameter  (millimeters) ,  and  low  power  (tens  of  milli¬ 
watts)  applications.  The  difficulty  in  the  use  of  stable  resonators  with 
large  apertures  and  high  power  will  become  apparent  in  the  discussions  in 
later  sections  of  this  report. 

The  second  type  of  resonator  has  become  extremely  useful  as  seen  for  the 
lasers  of  interest  in  this  text.  This  is  called  the  unstable  resonator  and 
was  originally  proposed  by  Slegman  [1] .  These  resonators  offer  many  advantages 
for  large  aperture,  high  power  CW  and  pulsed  applications. 

The  development  of  pulsed  lasers  leads  one  to  use  the  positive  branch 
unstable  confocal  resonator.  Mode  formation  time  is  not  of  interest  in  the 
CW  laser  case,  since  the  transverse  modes  will  sort  themselves  out  in  a  short 
time  compared  to  the  time  of  laser  operation.  However,  with  the  pulsed_laser , 
which  generates  power  over  a  relatively  short  interval  (0.01  to  50  x  10 
seconds)*,  the  situation  is  different.  In  this  case  it  became  apparent  that 
the  far-field  "focused"  energy  pattern  frequently  does  not  represent  the 
quasl-steady-state  distribution  as  calculated  from  the  transverse  modes  of 
the  CW  case. 

The  pioneer  paper  of  Zemskov,  et  al .  [2],  describes  the  experimental 
result  of  using  condition*.  ;.ly  stable  (plane  parallel)  versus  unstable 
resonators  in  a  short  pulse  copper  vapor  laser.  In  this  and  subsequent 
papers  [3,4]  it  was  determined  that,  aside  from  other  difficulties,  the  stable 
resonator  takes  an  unacceptably  large  number  of  resonator  round  trips  for  the 
radiation  to  sort  itself  out  into  the  transverse  modes  characteristic  of  a 
diffraction  limited  laser.  A  further  result  was  that  to  obtain  minimum 
beam  divergence  it  is  necessary  to  use  an  unstable  "telescopic"  resonator 
with  convergence  properties  based  upon  the  Fresnel  number  and  magnification 
of  the  resonator.  From  using  geometric  analysis  of  the  resonator  it  was 
discovered  and  experimentally  verified  that  a  very  large  magnification 
(M  •*  250)  is  required  to  Insure  near  diffraction  limited  performance  of  a 
laser  with  a  pulse  length  in  tens  of  nanoseconds. 

The  modeling  of  the  unstable  resonator,  described  by  Zemskov  as  well  as 
a  subsequent  refinement  by  Anan'ev  [5],  uses  this  geometric  approach  to 


analyze  transverse  node  formation.  In  the  Anan'ev  work,  the  angular  extent 
of  the  radiation  from  the  spontaneous  source  field  is  decreased  to  a  diffrac¬ 
tion  limited  output  by  successive  passes  between  convex  and  concave  resonator 
mirrors — a  transient  angular  "pinching  down"  of  laser  light  to  produce  a 
collimated  beam.  The  equation  he  derives  for  the  number  of  round  trips  re¬ 
quired  to  produce  a  quasi— stable  transverse  mode  in  a  confocal  unstable 
resonator  is  given  by 


n 


Jin  (M) 


(1.1) 


where  a  is  the  length  of  the  side  of  a  square  active  medium;  1,  the  laser 
wavelength;  -f  1 ,  the  secondary  mirror  focal  length;  and  K,  the  resonator 
magnification.  The  mode  formation  time  is  then  given  as 


where  L  is  the  resonator  mirror  separation,  and  a,  the  speed  of  light. 

This  basic  formula,  with  minor  modifications  and  adaptations  to  other 
resonators,  has  been  the  cornerstone  of  transverse  mode  formation  time  theory 
in  unstable  confocal  resonators  and  defines  the  minimum  number  of  round  trips 
to  produce  steady-state  mode  operation. 

This  geometrical  modeling  may  be  simply  Illustrated  in  Figure  1.  A 
typical  positive  branch  unstable  resonator  is  shown  with  a  concave  primary 
mirror,  m2,  of  half  dimension  a/2,  and  a  convex  secondary,  or 


Figure  1.  Geometry  of  transverse  mode  formation. 

feedback  mirror,  m.,  of  half  dimension  ax .  The  separation  of  -the  two  mirrors 
is  L.  This  analogy  is  valid  in  the  regime  of  moderate  magnification  (1.32  < 
M  <  3)  and  large  ratios  of  mirror  separation  to  beam  diameter.  After  each 
round  trip  the  angular  distribution  of  the  radiation  is  decreased  by  the 
resonator  magnification,  l.e.. 
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where  Che  resonator  magnification  is  given  by 


(1.8) 


with  f 2 ,  -  f1  being  the  respective  primary  and  secondary  mirror  focal  lengths, 
and 


l  -  2f 
n  n 


(1.9) 


A  criterion  for  the  minimum  far-fleld  spot  size  occurs  when  6^  becomes 
equal  to  a  diffraction  angle,  l.e. , 
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(1.10) 


where  A  is  the  laser  wavelength.  This,  combined  with  equation  (1.7),  gives 
the  approximate  number  of  round  trips  n",  required  for  the  unstable  resonator 
to  reach  the  diffraction  limit,  l.e. , 
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which  gives  the  same  approximate  number  of  round  trips  necessary  for  con¬ 
vergence  as  compared  to  the  Anan'ev  formula. 


B.  Stead-State  Transverse  Modes  in  Resonators 
1.  Stable  Resonators 

The  first  lasers  proposed  were  constructed  with  stable  resonator 
cavities  [6, 7, 8, 9]  and  as  a  result,  the  theory  of  stable  resonator  modes  has 
been  well  established.  An  excellent  example  of  this  development  is  given  in 
the  case  of  the  symmetric  confocal  resonator  [10,11].  The  empty  cavity 
transverse  modes  for  this  resonator  may  be  obtained  in  closed  form  as 


’  H»  [(ix)1*]  Ht  •*»  [-  (n) <x>  +  y!)] 


(1.14) 


where  U^Cx.y)  is  the  amplitude  eigenfunction  solution  of  the  transverse 

electromagnetic  (TEM)  modes  of  the  resonator;  H  and  H„  are  Hermite  polynomials 

®  * 

of  order  m  and  £:  L  and  X  have  been  defined  previously.  All  modes  are  power 
orthogonal  and  may  be  expressed  as  Hermite-Gauss  functions  for  resonators  with 
rectangular  symmetry.  Figure  2  is  a  schematic  of  the  symmetric  confocal 
resonator  with  a  few  of  the  lower  order  modes  illustrated  in  Figure  3. 

Stable  resonators  are  suitable  for  lower  power  applications;  Inherent 
difficulties  exist  in  their  high  power  (high  energy)  use.  Ideally,  we  desire 
to  have  all  the  energy  contained  within  the  lowest  loss  transverse  mode,  with 
the  mode  filling  the  cavity  active  media.  This  would  insure  beam  propagation 
to  the  far  field  with  maximum  fluence  and  minimum  angular  divergence  and  far 
field  spot  size.  Stable  resonators  do  not  have  these  characteristics.  In 
high  power  practical  applications  the  lowest-order  (usually  Gaussian)  mode 
is  small  compared  to  the  transverse  active  media  dimensions,  with  most  of  the 
energy  partitioned  into  higher  order  modes  with  increasingly  larger  divergence. 
As  an  example  [12]  the  lowest-order  mode  of  the  symmetric  confocal  resonator 
has  a  spot  size  2W  given  by 

2V-  2F  (8i,  d-15) 


LX  ** 

with  2Wq  =  2  —  .  F(gx ,  g2)is  a  number  which  depends  upon  resonator 

parameters  and  is  on  the  order  of  unity  or  slightly  larger  for  stable 
resonators.  For  a  typical  high  energy  C02  laser,  L  ■  200  cm  and  the  lowest 


order  mode  diameter,  2Wq,  Is  0.52  cm.  This  may  be  compared  to  a  desired  beam 

of  10  cm  else.  One  could  consider  using  stable  resonator  mirrors  having 
large  radii  of  curvature  to  fill  the  mode  volume  of  a  typical  high  energy 
laser  resonator.  However,  mechanical  stability  of  the  resonator  alignment 
becomes  a  limiting  practical  problem.  The  mirror  stability  for  nearly  flat 
mirrors  in  the  resonator  has  been  shown  by  Siegman  [12]  to  be 
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(1.16) 


i  determined  from  mirror  separation,  L,  and 
te  g  value  of  each  mirror  is  given  as 


(1.17) 

(1.18) 

determine  whether  or  not  a  reasonator  is  de- 
rdingly,  for  a  stable  resonator. 
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round  trips  of  a  planar  resonator  (Fabry- 
lce  to  a  steady-state  solution  has  been 
trge,  even  for  a  small  Fresnel  number  where 
s  expected  [6] . 


1  resonator  does  not  have  the  problems  noted 
ctive  coupling  results  in  large  mode  volume 
hat"  near-field  fluence  distribution  with 
tion  properties  in  the  far-field.  We  will 
ive  branch  confocal  unstable  resonator  [1,14] 
resonator  consists  of  a  small  convex 
ncave  ‘primary"  mirror  set  at  a  common  focus, 
is  determined  by  the  active  media  cross 
one try.  Figure  4  shows  the  resonator  and 
eometric  axis  in  the  center  of  the  two 
xis  in  the  corner  of  the  feedback  mirror 


this  resonator  are  predetermined  by  cavity 
s.  The  focal  length  (f x)  of  the  feedback 
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Figure  4.  Geometry  of  the  \ 
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.3  that  the  orthogonality  relations  over  the  cross  sections  are  of  the 
:orm 


//U  U  dxdy  «  6  (1.24) 

n  m  nm 

rather  than 

f/V  *U  dxdy  <*  <5  (1.25) 

n  m  1  nm 

.rhere  6  is  the  Kronecker  delta,  and  the  total  power  cannot  be  expressed 
nm 

as  the  sum  of  the  powers  carried  by  the  individual  modes  U  ,  U  . 


For  mirrors  having  rectangular  symmetry  as  illustrated  above,  the  eigen¬ 
functions  are  separable  in  x  and  y  and  are  in  the  form  Un(x)  Um(y) .  Figure  5 
shows  a  calculation  of  amplitude  and  phase  of  UQ(x)  across  a  single  transverse 
dimension  [17].  The  diffractive  rippling  of  the  near-field  as  seen  in  the  inten- 

i  i 2 

|U0(x)|  gives  the  two-dimensional  near-field  burn  pattern  of  Figure  6  as  a  charac¬ 
teristic  "plaid"  appearance  which  increases  in  fine  structure  with  increasing 
Fresnel  number .  Figure  7  shows  the  far-field  pattern  for  the  near-field  distri¬ 
bution  of  Figure  6. 

Higher  order  modes  of  rectangular  unstable  resonators  have  been  calculated 
by  numerical  techniques  [17,18,19]  and  analytic  techniques  involving  asymptot¬ 
ic  expansion  of  the  integral  equation  [20,21].  An  interesting  result  is  seen 
when  one  plots  the  power  losses  {l-|An|2}  for  an  unstable  resonator.  Figure  8 

[18]  shows  the  four  lowest  loss  modes  versus  equivalent  Fresnel  number. 
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(1.26) 


where  M  is  the  resonator  magnification  and  N^,  is  defined  as  the  tube  Fresnel 
number , 


(a/2)2 
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(1.27) 


It  is  interesting  to  note  that  at  low  equivalent  Fresnel  number  the  modes 
are  entwined  in  a  characteristic  oscillating  pattern  with  the  lowest  loss 
mode  (6)  eventually  separating  from  the  rest  at  higher  Fresnel  numbers. 

In  the  region  of  low  equivalent  Fresnel  number  (F  <  100) ,  the  most 
interesting  results  derived  from  detailed  theoretical  studies  [16]  are: 

(1)  the  lowest  loss  mode  is  always  symmetric; 

(2)  the  greatest  mode  loss  separation  occurs  at  F  "integer  +  .5; 
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Figure  6.  Near-field  burn  pattern,  M  ■  1.32,  NT  «*  279.  Vertical 

gaps  are  caused  by  the  feedback  mirror  holding  struts. 
The  pattern  is  actual  size. 
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(3)  equal  mode  loss  for  two  modes  occur  at  mode  crossing  points 
of  F  -  integer  +  .875.  Although  |Xn|  is  equal  for  the  two 

modes,  the  phase  is  distinct  and  different  for  each  mode; 

(4)  antisymmetric  modes  are  possible. 

The  lowest  loss  and  higher  loss  modes  may  be  obtained  by  a  numerical 
technique  using  Fourier  transform  methods  to  simulate  wavefront  propagation 
within  the  resonator  [22,23].  This  kind  of  calculation  utilizes  a  self- 
consistent  approach  in  which  the  wavefront  of  the  n+1  round  trip  resembles  the 
wavefront  of  the  nth  round  trip  to  an  arbitrarily  small  difference.  These 
techniques  have  been  modified  to  obtain  a  theoretical  basis  for  the  transverse 
■node  formation  work  described  in  this  text.  A  description  of  steady-state  mode 
calculations  in  unstable  resonators  by  the  Fourier  transform  approach  follows: 
the  lowest  order  mode  is  calculated  from  the  wave  equation,  which  is  a  dif¬ 
ferential  equation  analog  of  the  integral  eigenequation  (1.23)  with  the  appro¬ 
priate  boundary  conditions  of  the  unstable  resonator.  The  wave  equation  may 
be  obtained  directly  from  Maxwell's  equations  in  a  homogeneous  media  [24]: 


X  E  - 

3B 
‘  at 

(1.28) 

X  H  « 

j  +  M 

J  +  at 

(1.29) 

D  - 

e  E 

(1.30) 

B  - 

y  h 

(1.31) 

where  E  is  the  electric  field,  H  is  the  magnetic  field,  B  is  the  magnetic 

induction,  D  is  the  electric  displacement  and  J  is  the  current  density,  e  and 
y,  the  electric  permittivity  and  magnetic  permeability  of  the  medium, 
respectively. 

The  customary  paraxial  assumptions  of  light  propagation  in  a  homogeneous, 
isotropic  medium  are  made,  i.e. , 

7  (V  •  E)  -  0  (1.32) 
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with  and  -||  the  propagation  components  of  the  wave  equation  along  the 

Z-axis  and  k,  the  propagation  number.  It  has  been  shown  that  these  approxi¬ 
mations  are  valid  to  zero  order  derivation  of  Maxwell's  equations  for  in¬ 
homogeneous  media  with  focusing,  diffraction,  and  gain,  e.g. ,  active  lasting 
media  [25]. 


To  derive  the  wave  equation,  we  begin  with  the  vector  Identity 


(1.34) 


V  X  V  X  E  -  V  (V*E)  -VE  . 

Using  Equation  (1.32)  this  further  reduces  to 


V  X  V  X  E  -  -V2E 


Incorporating  Equations  (1.28)  and  (1.31)  into  (1.35): 


(1.35) 
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Now,  Equation  (1.29)  is  substituted  into  Equation  (1.36)  with  the  assumption 
that  no  electrical  currents  are  flowing  in  the  media,  i.e. ,  J  ■  0.  This  gives 
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Finally,  using  the  relation  of  Equation  (1.30),  we  have 


-y  92D  _  -ye  9^E 
at7  3t2 


(1.38) 


and  therefore,  since  ye 
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The  vector  equation  may  be  simplified  by  assuming  a  plane  polarized 
wave-front 


E  (x,y,z)  -  U(x,y , z)  E 


ilDt 


(1.40) 


Substituting  E  into  Equation  (1.39)  removes  the  temporal  and  vector  dependence 
to  give  the  Helmholtz  Equation 


vz  u  -  -  ~  u 


(1.41) 


with 


y.2  -  U 
c2 


(1.42) 


Considering  the  paraxial  movement  of  light  through  the  medium,  U 
may  be  decomposed  into  a  slowly  varying  component  multiplied  times  the 
characteristic  plane  wave  propagation  at  a  single  frequency: 
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U(x,y,z)  -  \p(x,y,z)  e1 


(1. 


Substituting  this  assumption  for  U  into  Equation  (1.41)  we  obtain 


and  applying  the  paraxial  approximation 


(1. 


|!i 

«  k 

3z* 

3z 

(1.: 


we  obtain  the  wave  equation  of  propagation 


0 


(!•' 


The  Fourier  transform  lends  itself  readily  to  the  solution  of  Equation 
(1.45).  To  se  this,  we  examine  a  basic  property  of  these  transforms.  The 
Fourier  transform  of  ip  in  the  x  dimension  is  given  as 


T 


x 


(W 


e  *px  \p(x)  dx 


(l.< 


then. 


’„(£)  -r  «-ip*  i  - 


a.' 


Integrating  by  parts  we  get 


*48)  ■  I  ”  e'lp*  * +  ip  /”  ''lp,*“  • 

Since  ^  (x,y,z)  vanishes  at  infinity  we  get 
Tx  S  "  1  p  TxW  "  * 


(!.' 


(1.- 


by  a  similar  procedure 


1 


Taking  the  Fourier  transform  of  x  and  y  in  Equation  (1.45)  we  have 

-(p2  +  q2)  +  2  ik  -  0  .  (1.51) 

The  transformed  wavefront  ip  may  be  readily  propagated  from  zq  to  z^  by 
solution  of  Equation  (1.51): 
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(p2  ±  q2) 

2  lk 


dz 


(1.52) 


in  ip  (p.q.z^)  -  in  (p,q,zo) 


2  .  2 

p  +  q 

2  lk 


(z 


-  zo> 


Therefore,  by  the  property  of  logarithms, 


(1.53) 


M.q.z  ) 


-i(p2  + 


^  (p.q.z  ) 
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The  propagated  complex  amplitude  may  then  be  readily  obtained  as 

»\i 

(x,y,z  )  -  T  _1[^(p,q,z  ) ]  .  ± 


(1.54) 


1.55) 


Thus,  an  arbitrary  wavefront  may  be  launched  within  the  confines  of  an  unstable 
resonator  by  using  this  method.  The  direct  and  inverse  Fourier  transforms  are 
calculated  numerically  as  is  described  in  a  later  chapter. 

Far-fleld  calculations  and  experiments  show  that  the  Fraunhofer  diffrac¬ 
tion  pattern  from  the  unstable  resonator  is  essentially  the  far-field  pattern 
obtained  from  a  constant  phase-constant  amplitude  wavefront  illuminating  the 
geometric  aperture  [21,26,27,28].  The  far-fleld  pattern  of  this  case  has  been 
calculated  for  centered  as  well  as  completely  decentered  obscurations  for  the 
empty  cavity  case.  Calculations  using  gain  appear  to  give  similar  results 
[29,30]. 


C.  Discussion  of  Relevance  of  the  Problem 

The  goal  of  this  research  was  to  examine  the  temporal  evolution 
of  transverse  modes  in  moderate  magnification  (1.32  <.  M  <.  2.70),  positive 
branch  confocal  unstable  high  energy  laser  resonators.  Toward  this  end,  the 
far-fleld  irradiance  (fluence)  distribution  has  been  observed  in  arbitrary 
time  slices  during  and  following  the  period  of  mode  formation.  The  geometrical 
theory  of  Anan'ev  was  compared  to  experimental  results  and  statistical  analyti¬ 
cal  modeling  in  this  work. 

The  phenomenon  of  mode  formation  is  illustrated  schematically  in  Figure  9. 

A  resonator  cevlty  is  shown  with  active  media  on  the  left  with  a  resultant 
far-field  irradiance  distribution  on  the  right.  At  threshold  (t  -  o),  radiation 
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Figure  9.  Mode  formation  observed  in  the  far-field. 

is  emitted  spontaneously  in  all  directions,  the  far-field  lrradiance  being 
exemplified  by  low  level  noise.  At  a  later  time  during  mode  formation 
(t  -  200  nsec),  a  far-field  pattern  is  beginning  to  form,  represented  by 
skewed  and  displaced  peaks.  Finally,  at  a  time  longer  than  that  of  mode 
formation  (t  -  1000  nsec) ,  a  far-field  lrradiance  distribution  has  stabilized 
to  form  the  distribution  characteristic  of  steady-state  mode  generation. 

A  situation  can  readily  exist  in  which  the  spontaneous  emission  and  gain 
in  the  resonator  are  high  at  threshold.  In  this  case  a  large  initial  intensity 
builds  up  and  saturation  of  flux  can  occur  during  the  period  of  predicted  mode 
formation.  If  this  is  Indeed  the  case,  a  sizeable  percentage  of  the  total 
energy  output  from  a  pulsed  laser  may  be  lost  by  mode  formation  effects  from 
the  central  lobe  of  far-field  pattern. 

Notable  high  energy  pulsed  lasers  which  have  short  start-up  times  and  fall 
under  the  category  described  above,  are  rare  gas  halide  (excimer)  lasers, 
metal  vapor  lasers,  and  chemical  lasers.  Some  work  has  been  done  on  mode 
formation  and  beam  quality  for  copper  vapor  and  excimer  lasers  utilizing  un¬ 
stable  resonators  [2,4,31].  These  experiments  examined  resonators  with 
magnifications  which  were  relatively  large  (>.  8)  and  thus  had  mode  formation 
times  of  the  order  of  tens  of  nanoseconds.  Several  authors  have  derived 

18 


* 


i 


estimates  of  mode  formation  time  based  upon  the  geometric  theory  approach  [2, 
4,5,31,32],  but  experimental  data  is  lacking  in  the  lovermagnif ication  region 
of  interest  for  high  energy  excimer  and  chemical  lasers  (1  <  M  <  3) . 

Resonators  in  this  regime  are  predicted  to  have  mode  formation  times  in  the 
hundreds  of  nanoseconds. 

The  prospect  of  mode  formation  times  of  a  few  hundred  nanoseconds  does 
not  bode  well  for  high  energy  excimer  and  chemical  lasers.  Near-field  power 
turn-on  in  excimer  systems  has  been  observed  in  the  50  to  200  nanosecond 
regime  [33].  This  has  been  confirmed  independently  by  the  author  and 
J.  Oldenettel  in  unpublished  work.  Pulsed  chemical  lasers  are  seen  to  have 
10  to  40  lines,  which  turn  on  during  a  3  microsecond  pulse;  each  line  being 
active  for  one  to  a  few  hundred  nanoseconds  [34].  Thus,  mode  formation  can 
create  a  serious  problem  for  these  lasers  with  delivery  of  a  divergent  beam 
to  the  far-field  in  the  initial  portion  of  the  pulse. 

The  problem  may  be  succinctly  stated  as  follows:  certain  common  lasers 
used  in  pulsed  high  energy  laser  systems  may  be  inefficient  in  delivery  of 
far-field  energy  because  of  finite  mode  formation  time.  A  confirmation  of 
the  geometrically  predicted  limit  has  been  attempted  in  this  treatise  through 
experimental  investigation  and  statistical  computational  modeling. 

II.  EXPERIMENTAL  RESULTS 

A.  Description  of  Lasers  Used  in  This  Work 
1.  Types  of  Lasers 

It  is  desirable  to  validate  transverse  mode  formation  time  de¬ 
fined  by  Equation  (1.2)  over  the  range  of  values  of  variables  encountered  in 
practical  high  energy  lasers.  As  discussed  in  Section  1,  this  phenomenon  may 
be  of  critical  importance  in  the  design  and  implementation  of  a  high  energy 
laser  system.  To  experimentally  measure  mode  formation  an  ideal  laser  would 
have  the  following  characteristics: 

a.  high  initial  spontaneous  emission 

b.  high  gain 

c.  good  active  medium  quality 

d.  variable  magnification 

e.  variable  Fresnel  number 

Characteristics  a,  b,  and  c  refer  to  the  lasing  medium.  It  would  be  desirable 
for  the  laser  to  "turn  on"  quickly  after  threshold  and  saturate  so  that  near¬ 
field  output  power  is  nominally  constant  and  far-field  measurements  reflect 
mode  formation  properties.  In  addition,  the  active  medium  must  not  contribute 
time  varying  phase  distortions  which  would  confuse  far-field  measuresients. 

Parts  d  and  e  relate  directly  to  Equation  (1.2)  through  the  resonator  geometry. 
In  regions  of  low  magnification,  M  becomes  the  predominant  factor  in  the 
denominator.  Fresnel  number  may  vary  over  a  large  range  but  is  not  as  important 
as  M  because  of  the  logarithmic  relationship. 
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Figure  10.  Pulsed  high  energy  laser  schematic. 

The  COj  laser  consisted  of  three  basic  units:  the  electron  beam  (e-beam) 
supply  which  ionized  the  laser  gas,  the  resonator  cavity  containing  the 
active  laser  media  with  unstable  confocal  resonator  optics,  and  the  discharge 
supply — a  DC  source  which  maintained  a  constant  voltage  across  the  laser  cavity 
and  determines  the  gain  in  the  active  lasing  medium.  The  e-beam  was  a  cold 
cathode  emitter  type;  electrons  were  emitted  from  the  edges  of  .001  in.  thick 
tantalum  foil  strips  on  the  cathode  assembly.  The  e-beam  pulse  length  was 
variably  controlled  up  to  4  microseconds  with  a  nominal  voltage  of  180  kV. 

The  e-beam  gun  chamber  was  held  at  low  pressure  (y  1  x  10~6  Torr)  and  was 
separated  from  the  laser  cavity  by  a  .001  in.  thickness  titanium  foil.  The 
current  density  of  the  e-beam  on  the  cavity  side  of  the  foil  was  approximately 
0.5  amperes/cnr  for  a  beam  dimension  of  10  x  100  centimeters.  The  laser  cavity 
had  an  active  volume  of  10  x  10  x  100  centimeters  separated  on  the  longitudinal 
axis  by  unstable  confocal  positive  branch  resonator  mirrors.  For  these  experi¬ 
ments,  the  resonator  mirror  pair  had  a  magnification,  M,  of  1.32  at  a  confocal 
spacing  of  289  cm.  The  convex  mirror  (feedback  mirror),  which  is  7.5  cm  in 
diameter,  gave  an  output  beam  of  9.9  cm  through  a  NaCl  output  window.  The 
concave  reflector  (primary  mirror)  was  oversized  at  16  cm  diameter.  The 
mirrors  were  separated  by  L  -  337.5  cm  to  give  focusing  at  136  meters,  the 
nominal  position  for  far-fleld  analysis  of  the  laser  beam.  In  terms  of 
resonator  g  parameters  defined  by  Equation  (1.17),  focusing  in  terms  of  units 
of  mirror  separation,  L,  as  per  Siegman  [12]  is  given  as 


T-l)  + 


8^)  8x 


(2.1) 


r  ,#■ 


TABLE  I.  LASER  SYSTEM  PARAMETERS 


Parameter 

CO2  Laser 

XeF  Laser 

Laser  Wavelength 

10. 6u 

0.3532u 

Pulse  Length  (In  Experiments) 

0.02-1.0  y  sec 

0.2-0. 8  ysec 

Pulse  Energy  (In  Experiments) 

30  J 

35  J 

Laser  Gas  Mix  (In  Experiments) 

1:1  C02:He 

1:3:1140  NF,:X  :N 
3  e  i 

Cavity  Pressure  (In  Experiments) 

456  Torr 

1520-3040  Torr 

Active  Media  Dimensions 

10x10x100  cm 

10x20x200  cm 

Resonator  Magnification 

1.32 

1.32,  2.7 

Mirror  Separation 

338  cm 

262  cm 

Output  Beam  Size 

10x10  cm 

10x20  cm 

Geometric  Out coupling 

0.426 

0.426,  0.863 

Beam  Shape  (Obscuration) 

Decentered 

Centered 

Tube  Fresnel  Number 

279 

2700;  10,800 

Mode  Formation  Time 

489  nsec 

513,  168  nsec 

E-Beam  Voltage 

180  keV 

300-500  keV 

Discharge  Voltage 

0-50  keV 

0 

Guide  Fielc  Current 

0 

300  A 

where  in  this  case  LF'  ■  136  meters.  The  mirror  radii  of  curvature  may  be 
determined  by  Equations  (1.8),  (1.9),  and  (1.21)  and  used  in  Equation  (1.17). 
The  laser  medium  was  a  mixture  of  equal  parts  of  carbon  dioxide  and  helium  at 
456  Torr  pressure.  Transverse  pumping  of  the  laser  gas  was  performed  by  a  DC 
discharge  carbon/graphite  anode  of  somewhat  larger  dimensions  than  the 
10  x  100  cm  e-beam  area  facing  the  e-beam  cathode  foil  with  11  cm  separation. 
The  anode  was  preset  in  a  range  of  40  -  50  kV  during  the  e-beam  pulse.  The 
e-beam  determines  the  approximate  laser  pulse  length. 

The  discharge  supply  provided  the  basic  pump  mechanism  of  the  CO 2  laser . 

This  is  a  source  of  DC  current  which  is  arbitrarily  chosen  for  a  given  pulse. 

A  DC  supply  charges  a  set  of  capacitors  which  are  then  discharged  ("crowbarred") 
after  cessation  of  the  e-beam  pulse.  The  actual  C0«  laser  assembly  used  for 
these  measurements  is  shown  in  Figure  11,  without  tne  resonator  mirrors. 
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Figure  11.  CO.  laser. 


The  XeF  laser  was  very  similar  to  the  CO 2  laser  in  construction.  However, 

the  e-beam  was  a  higher  voltage  (300  -  500  kV)  and  the  dimensions  were  larger 
(20  cm  high  by  200  cm  long)  than  those  of  the  CO 2  laser.  The  basic  pumping 

mechanism  in  the  XeF  laser  was  through  direct  e-beam  initiation;  no  DC  dis¬ 
charge  was  used.  The  mirrors  were  set  at  the  confocal  spacing  of  262  cm.  Two 
resonator  mirror-pairs  were  chosen  for  the  experiments:  an  M  ■  1.32  set  and 
M  =  2.7  set,  which  gave  a  reasonable  range  of  mode  formation  times  according 
to  Equation  (1.2).  The  cavity  pressures  in  this  laser  were  higher  than  in  the 
C02  laser,  being  in  the  range  of  2  -  4  atmospheres.  The  typical  active  lasing 

medium  was  10  x  20  x  200  cm  in  volume  with  a  gas  mixture  ratio  of  1  to  3  to 
1140  for  nitrogen  tribluoride  (NF3)  to  xenon  to  neon.  A  magnetic  guide  field 

was  produced  by  a  pair  of  Helmholtz  coils  surrounding  the  cavity  for  e-beam 
confinement.  The  laser  used  for  these  measurements  is  shown  in  Figure  12. 

2.  Pumping  Mechanisms 

The  value  of  small  signal  gain,  y  ,  is  determined  by  the 
following  equation  [ 39 ] 


Y 
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A2g(v) 


r|28nt 


spont 


(2.2) 


where  X  is  the  wavelength  of  the  laser  line,  g(v)  is  the  absorption  normalized 
lineshape,  n  is  the  index  of  refraction,  tSpont  the  spontaneous  emission 

lifetime  of  the  medium  N2  and  Nx  are  the  upper  and  lower  laser  level  populations; 
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Figure  13.  Pumping  scheme  for  the  pulsed  C02  laser  [43] 

(Nominal  3:2:1,  He :N  :C0  Mix). 

2  2 

The  lower  level  (100)  is  quickly  deactivated  by  collisional/radiative 
processes  to  the  ground  state.  Collisional  deactivation  is  increased  with 
:he  aid  of  the  He  diluent. 

The  gain  in  the  transversely  pumped  high  pressure  C02  laser  has  been 

3hown  to  be  proportional  to  E/N  [44,45,46]  where  E  is  the  transverse  discharge 
Field  (in  volts/cm)  and  N  is  the  number  of  molecules  per  cubic  centimeter 
(proportional  to  pressure) . 

The  XeF  laser  kinetics  are  not  straightforward  as  compared  to  the  C02 

Laser.  A  number  of  complicated  paths  are  possible  for  activation  of  the  upper 
Laser  level  represented  by  XeF*.  These  are  shown  in  Figure  14  [47].  The 
najor  energy  flow  paths  in  the  Figure  are  given  as: 

A 

a.  pumping  electrons  Ne  -*■  NeF  ■>  F  ■>  Xe  -*■  XeF 

b.  pumping  electrons  -*■  Ne  -*•  Xe  -+•  XeF 

c.  pumping  electrons  -*•  Xe+  •>  XeF* 

A  'fc 

d.  pumping  electrons  -*■  Xe  -*■  XeF 

rhese  expressions  are  not  equations,  but  are  shown  only  as  examples  of 
energy  flov  to  illustrate  the  complexity  of  XeF  pumping.  With  the  pressures 
and  diluent  ratios  of  neon  or  argon  described  in  this  work,  the  B  -*■  X  state 
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(3532  A)  is  pumped  [48].  The  lower  level  of  this  state  was  originally 
thought  to  be  purely  dissociative.,  i.e. , 


2XeF*  -*•  2Xe  +  F2  +  v 


(2.5) 


However,  a  bound  state  of  XeF  exists  in  the  lower  level  below  approximately 
500 °C,  which  reduces  the  overall  population  inversion  and  gain. 

B.  Small  Signal  Gain  Measurements  in  C02  Lasers 


The  time  of  transverse  mode  formation.  Equation  (1.2),  is  measured 
from  the  time  at  which  threshold  gain,  g  ,  is  achieved,  i.e. ,  feedback  power 
over  one  round  trip  pass  equals  out  coup?ed  power.  The  method  for  measuring 
threshold  gain  is  illustrated  by  Figure  15.  This  is  shown  for  C02;  similar 

measurements  are  possible  for  XeF  with  substitute  of  a  suitable  detector  and 
probe  laser  and  will  be  explained  later.  The  experiment  involves  a  single  pass 
C02  probe  laser  turned  to  the  high  energy  laser  line  and  a  high  speed  detector 

to  make  small  signal  gain  measurements.  DC  levels  of  I  *  0  and  I  are  deter¬ 
mined  and  the  high  energy  laser  is  pumped  under  experimental  operating 
conditions.  The  temporal  history  of  the  power  increase,  dl/dt,  is  measured  as 
a  function  of  time  for  each  pumping  condition.  The  expression  for  gain  in  a 
single  pass  is  given  as 


HO  -  eg0(0* 


(2.6) 


where  £  is  the  active  medium  length.  In  actual  high  energy  laser  operation  the 
threshold  for  a  single  round  trip  in  the  resonator  occurs  at 


-e  2*o(t)* 


(2.7) 


where  l'(t)  and  I  '  are  intensities  measured  under  saturable  gain  conditions, 
o 

Therefore , 
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and  threshold  gain  occurs  when  l/Io  is  equal  to  the  resonator  magnification. 

The  increase  in  I  is  related  to  a  reference  signal  which  is  available  during 
high  energy  lasing.  This  reference  signal  is  the  current  produced  by  the 
electron  beam  entering  the  laser  cavity  (e-beam  postfoil  current) . 

Figure  15  shows  the  mechanics  of  making  gain  measurements  with  the 
Systems,  Science,  and  Software  (S3)  pulsed  CO  laser  with  optics  removed.  A 

CW  CO  laser  set  on  the  pulsed  laser  line  (P-20)  is  transmitted  through  the  S3 

2 

laser  to  an  RF  shielded  fast  response  mercury-cadmium-tellur ide  detector.  The 
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Figure  15.  Gain  measurement  schematic. 

output  signal  is  displayed  on  a  7844  dual  beam  oscilloscope.  A  spectrograph 
checks  the  P-20  line  immediately  before  and  after  S3  firing  with  the  beam 
deflected  with  a  removable  mirror.  A  helium  neon  laser  is  used  for  visible 
alignment  of  all  components,  the  helium  neon  beam  following  the  path  of  the 
CW  CO 2  beam.  Zero  level  intensity  (I  =  0)  is  obtained  by  closing  a  shutter 
in  front  of  the  laser  and  allowing  the  oscilloscope  to  sweep  immediately 
before  the  test. 

Figures  16  through  20  show  small  signal  gain  measurements  as  a  function 
of  time  for  various  gain  (discharge  voltage)  settings  at  456  Torr  cavity 
pressure  and  100  cm  gain  length.  The  bottom  trace  is  the  e-beam  postfoil 
current,  the  trace  one  division  above  is  the  1=0  line,  and  the  upper 
trace  is  I  *  I(t)  with  the  flat  portion  to  the  left  giving  I  =  IQ.  Threshold 

gain  is  measured  from  the  rise  of  postfoil  current  ot  the  point  at  which 
I/I  =1.32.  Figure  20  is  interesting  in  that  I(t)  is  shown  over  the  nominal 
duration  of  the  e-beam  pulse  in  which  I/I  exceeds  40  (small  signal  gain  > 
3.7%/cm.) 


Figure  16.  Small  signal  gain  at  40.8  kV  discharge  voltage. 


Figure  17.  Small  signal  gain  at  46.9  kV  discharge  voltage 


Figure  18.  Small  signal  gain  at  48.9  kV  discharge  voltage 


Figure  19.  Small  signal  gain  at  50  kV  discharge  voltage 


Figure  21.  Laser  output  et  40.8  kV  discharge  voltage. 


Figure  22.  Laser  output  at  46.8  kV  discharge  voltage 


Figure  23.  Laser  output  at  48.8  kV  discharge  voltage 


Figure  20.  Small  signal  gain  at  47.9  kV  discharge  voltage  (1q  ■  8mV) . 

Figures  21  through  24  show  high  power  output  from  the  S3  laser  with 
unstable  resonator  optics  installed.  The  oscilloscope  traces  were  taken  with 
a  7441  photon  drag  detector  in  the  near-field  in  the  upper  portion  of  the  beam 
midway  between  anode  and  cathode.  Again,  the  postfoil  current  is  being  used 
as  a  reference  pulse  and  is  the  lower  trace.  The  upper  trace  is  the  photon 
drag  detector  response.  These  photographs,  along  with  Figures  16  through  19 
of  small  signal  gain,  determine  the  time  of  threshold  relative  to  the  onset  of 
high  peak  power  from  the  laser  by  using  the  e-beam  postfoil  current  as  a 
reference  marker.  Increasing  gain  shortens  the  time  between  threshold  and 
high  peak  output  power  from  the  S3  laser.  Transverse  mode  formation  is  mea¬ 
sured  from  the  threshold  to  a  point  required  to  produce  a  far-field  burn 
pattern  in  these  experiments.  This  occurs  at  a  time  determined  by  the  onset 
of  high  peak  power.  The  photon  drag  detector  samples  2.3%  of  the  total  power 
in  Figures  21  through  24.  In  the  far-field  measurements  of  mode  formation 
time  to  be  presented  later,  the  detector  sampled  approximately  3%  of  the  total 
power.  The  determination  of  burn  pattern  formation  time  from  these  figures 
assumes  a  rise  in  power  of  40  mV  on  the  scale  of  Figures  21  through  24, 
referenced  to  e-beam  postfoil  current  rise  time.  Table  II  summarizes  the 
results  of  burn  pattern  formation  time  and  threshold  gain  time  as  a  function 
of  discharge  voltage.  The  last  column  of  numbers  represents  the  actual  mode 
formation  time  to  obtain  a  far-field  pattern  burn  and  is  the  Anan'ev  time  for 
the  following  CO  experiments. 

C.  Far-Field  CO z  Burn  Patterns 

The  experimental  arrangement  for  performing  transverse  mode  forma¬ 
tion  is  shown  in  Figure  25  and  was  placed  at  136  meters  from  the  S3  laser. 

The  laser  beam  focus  is  nominally  at  136  meters.  The  beam  is  propagated  an 
additional  distance  of  2.3  meters  to  a  zinc  selenide  beamsplitter.  Immediately 
behind  the  beamsplitter  is  a  plastic  sheet  attenuator  (3M  infrared  trans¬ 
parency)  of  <  30%  transmission  and  Polaroid  Sheet  Px  for  recording  long  pulse 
(1.2  y  sec)  burn  patterns.  This  sheet  is  unexposed,  developed  Polaroid  type 
52  film.  A  front  surface  reflection  from  the  beamsplitter  is  propagated  3.3 
meters  to  concave  mirror  mz  which  focuses  the  beam  at  1.2  meters  with  air 
breakdown  at  this  position,  P2.  The  air  breakdown  is  used  as  a  fast  plasma 
switch  to  clip  all  but  the  leading  edge  of  the  pulse.  This  is  important  in 
attempting  to  obtain  a  burn  pattern  at  P4  in  the  regime  of  mode  formation  as 


Figure  24.  Laser  output  at  50  kV  discharge  voltage. 

TABLE  II.  FAR -FIELD  BURN  PATTERN  FORMATION  TIME  VERSUS  DISCHARGE 


VOLTAGE. 

DISCHARGE 
VOLTAGE  (kV) 

BURN  PATTERN 
FORMATION  TIME  (ns) 

THRESHOLD 

GAIN  TIME  (ns) 

MODE  FORMATION 
TIME  (ns) 

40.8 

740 

220 

520 

46.8 

650 

190 

460 

48.8 

640 

180 

460 

50.0 

620 

180 

440 

Figure  25.  Experimental  schematic  to  obtain  mode  formation  data 


predicted  by  Anan'ev.  A  long  pulse  burn  pattern  would  have  integrated 
energy  long  after  the  transverse  mode  was  expected  to  have  been  formed.  After 
air  breakdown  at  P2  occurs,  the  beam  propagates  an  additional  2.2  meters  to 
concave  mirror  m2  where  the  geometrical  outline  of  the  laser  cavity  is  re¬ 
imaged  at  P3,  a  distance  of  1.8  meters  from  m2 .  At  this  position  an  aperture 
stop  is  inserted.  This  is  an  oversized  L-shaped  hole  which  contains  the  image 
of  the  near-field  laser  aperture.  The  purpose  of  the  aperture  stop  is  to 
block  rays  of  light  transmitted  after  air  breakdown  at  P2.  A  beamsplitter  is 
positioned  2  cm  downstream  of  the  spatial  filter.  Ninety  percent  of  the 
clipped  pulse  forms  a  burn  pattern  at  P4.  Ten  percent  of  the  beam  is  inter¬ 
cepted  by  a  photon  drag  detector  at  P5.  This  signal  is  received  by  an  os¬ 
cilloscope  which  gives  the  temporal  shape  of  the  pulse. 

The  objective  of  these  experiments  was  to  obtain  results  as  early  as 
possible  in  the  temporal  formation  of  a  transverse  mode.  Toward  this  end 
the  E/N  ratio  of  the  1:1  C02:He  laser  was  maximized  by  setting  the  gas  mixture 
to  456  Torr.  N2  was  not  used  because  of  its  long  N2-C02  transfer  time.  The 

discharge  voltage  was  maximized  in  the  regime  of  40  to  50  kV  to  provide  for 
high  initial  gain;  the  e-beam  pulse  width  was  minimized  to  protect  against 
laser  device  damage;  and  the  air  breakdown  plasma  switch  was  used  to  produce 
burn  patterns  representative  of  an  early  time  in  the  pulse.  The  time  of  burn 
pattern  formation  from  threshold  as  a  function  of  discharge  voltage  is  given 
in  Table  II. 

Three  sets  of  figures  are  representative  of  the  burn  patterns  taken  at  P4 
of  the  clipped  pulse  in  the  far-field.  Set  I  (Figures  26  through  29)  was 
taken  at  a  discharge  voltage  of  40.8  kV  which  gives  a  burn  pattern  formation 
of  520  nanoseconds  -  somewhat  outside  the  nominal  Anan'ev  mode  formation 
prediction  of  489  nanoseconds.  Set  II  (Figures  30  through  33)  is  burns  taken 
with  a  discharge  voltage  of  47.8  kV,  representing  a  burn  pattern  formed  at 
460  nanoseconds.  Set  III  (Figures  34  through  37)  was  taken  at  the  voltage 
limit  of  the  tests,  with  a  discharge  voltage  of  kV  giving  the  shortest  burn 
pattern  formation  time  of  440  nanoseconds. 

An  inspection  of  the  three  sets  of  figures  shows  a  successive  blurring 
and  subsequent  loss  of  resolution  in  the  far-field  pattern — a  result  attri¬ 
butable  to  mode  formation.  Set  III,  in  face,  has  been  blurred  to  an  extreme 
extent.  Only  the  tip  of  the  central  maximum  has  sufficient  fluence  to  burn  the 
Polaroid.  This  interpretation  was  corroborated  by  the  photon  drag  detector 
signal  at  P3<  Three  representative  temporal  oscilloscope  photographs  are 
shown  in  Figures  38  through  40  for  40.8  kV,  47.8  kV,  and  50  kV  discharge 
voltages,  respectively.  The  total  energy  obtained  from  integration  of  the 
photon  drag  temporal  response  in  each  figure  is  approximately  the  same.  One 
can  conclude  that  the  energy  was  spread  (blurred)  over  a  larger  area,  as 
anticipated  from  the  dynamics  of  mode  formation.  However,  the  energy  was 
collected  in  a  sufficiently  large  area  by  the  photon  drag  detector  to  compen¬ 
sate  for  the  blurring  of  the  far-field  pattern,  and  recorded  the  same  energy 
for  Figures  38  through  40. 

Figure  41  is  representative  of  all  long  pulse  burn  patterns  taken  at  Pj^ 
simultaneously  with  the  patterns  taken  at  P4.  The  purpose  of  this  burn  is  to 
Illustrate  the  quasi-steady-state  far-field  pattern  from  the  laser.  The  burn 
pattern  in  this  figure  has  formed  at  least  1,000  nanoseconds  after  threshold. 
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Figure  26.  Far-field  pattern  taken  at  40.8  kV  discharge 
(Shot  number  7174)  7X  enlargement.) 


Figure  28.  Far-field  pattern  taken  at  40.8  kV  discharge 
(Shot  number  7207)  7X  enlargement. 


Figure  29.  Far-field  pattern  taken  at  40.8  kV  discharge 
(Shot  number  7154)  7X  enlargement. 
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Figure  30.  Far-field  pattern  taken  at  47.8  kV  discharge 
(Shot  number  7197)  7X  enlargement . 


Figure  31.  Far-field  pattern  taken  at  47.8  kV  discharge 
(Shot  number  7210)  7X  enlargement. 


Figure  32.  Far-field  pattern  taken  at  47.8  kV  discharge 
(Shot  number  7208)  7X  enlargement. 


Figure  33.  Far-field  pattern  taken  at  47.8  kV  discharge 
(Shot  number  7196)  7X  enlargement. 

i 
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Figure  34.  Far-field  pattern  taken  at  50  kV  discharge 
(Shot  number  7198)  7X  enlargement. 


Figure  35.  Far-field  pattern  taken  at  50  kV  discharge 
(Shot  number  7219)  7X  enlargement. 


37 


Figure  36.  Far-field  pattern  taken  at  50  kV  discharge 
(Shot  number  7202)  7X  enlargement. 


Figure  37.  Far-field  pattern  taken  at  50  kV  discharge 
(Shot  number  7201)  7X  enlargement. 


Figure  38.  Temporal  response  of  clipped  pulse  at  40.8  kV 
discharge  (Shot  number  7207). 


Figure  39.  Temporal  response  of  clipped  pulse  at  47.8  kV 
discharge  (Shot  number  7197) . 


Figure  40.  Temporal  response  of  clipped  pulse  at  50  kV 
discharge  (Shot  number  7201). 


Figure  41.  Long  pulse  burn  pattern  from  laser  at  50  kV 
discharge  (Shot  number  7201) . 

The  burn  pattern  at  this  position,  except  for  overall  fluence,  is  independent 
of  discharge  voltage  occurring  well  after  the  mode  formation  time  predicted 
by  Anan'ev  [ 5] . 

The  temporal  response  of  the  long  pulse  burn  compared  to  the  e-beam 
pulse  is  shown  in  Figure  42  for  a  discharge  voltage  of  47.8  kV.  The  long  pulse 
trace  has  been  electronically  shifted  back  600  ns.  as  compared  to  the  e-beam 
initiation  in  order  to  observe  the  characteristics  of  the  total  pulse. 


Figure  42.  Long  pulse  burn  temporal  response  (lower  trace) 
and  E-beam  current  (upper  trace)  . 

In  order  that  a  comparison  may  be  made  to  the  long  pulse  burn  of  Figure  41, 
an  ideal  steady-state  calculation  of  relative  fluence  (joules/cm2)  contours 
has  been  made,  and  is  shown  in  Figure  43.  This  is  a  good  approximation  to  the 
Fraunhofer  pattern  at  the  focus  of  the  laser  at  136  meters  and  is  the  pattern 
resulting  from  a  constant  phase,  constant  amplitude  "top  hat"  distribution 
across  the  laser  geometrical  aperture.  The  L-shaped  beam  from  the  laser  for 
these  experiments — a  result  of  complete  decentering  of  the  feedback  mirror 
(Figure  4) — gives  a  far-field  burn  (Figure  41)  which  is  notable  different  in 
shape  compared  to  the  centered  burn  pattern  (Figure  7)  for  this  resonator.  The 
decentered  configuration  was  chosen  for  two  reasons:  the  increase  of  Fresnel 
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Figure  43.  Ideal  far-field  pattern  at  136  meters. 


number  to  increase  mode  formation  time  and  the  ease  of  preparation  of  a  spatial 
filtering  mask,  which  for  the  decentered  resonator  is  an  L-shaped  aperture. 

The  relative  fluence  contours  in  Figure  43  are  adjusted  logarithmically  at 
five  per  decade  (1.58  X  per  contour).  This  is  done  to  give  a  realistic  com¬ 
parison  to  the  nonlinear  burns  of  Figures  26  through  37  and  Figure  41. 

Calorimetry  calibrations  were  attempted  at  Pj  and  P5  and  in  the  laser 
near-field  using  the  Scientech  364  Joulemeter  (see  Figure  25) .  This  meter 
measures  the  temperature  rise  of  a  low  reflectance  thermopile,  causing  a 
change  of  voltage  output  that  is  interpreted  as  energy.  The  energy  collected 
at  Pi,  was  0.1  Joules — typical  number  for  the  burn  patterns  fo  Figures  26  through 
37,  which  have  been  enlarged  seven  times  in  the  photographs.  The  energy  was 
collected  over  twice  the  burn  dimensions  as  shown  on  the  photographs. 


By  comparison,  a  similar  collection  of  energy  without  the  plastic  atten¬ 
uator  at  P x  was  approximately  26  Joules  .  With  the  plastic  attenuator  in  place 
this  would  have  been  6-8  Joules  .  Figure  41,  which  was  taken  at  Px  has  about 
the  same  depth  of  burn  as  Figure  26,  which  was  taken  at  P4,  and  resulted  in 
the  same  recorded  fluence.  Taking  into  account  the  15%  reflection  per  side  of 
the  zinc  selenide  beamsplitter,  the  10:1  beamsplitter  near  P4,  and  the  7X 
demagnification  of  the  burn  at  P4;  P4  should  have  an  energy  reading  based  on  8 


joules  at  Pj 


of :  Joules  equals  0.1  Joules,  which  was  measured  at 


P4.  The  Type  52  Polaroid  sheets  used  to  obtain  the  burn  patterns  in  Figures 
26  through  37  and  Figure  41  had  a  minimum  sensitivity  of  approximately  0.15 
Joules/cm"  fluence.  The  fluence  for  Figures  26  through  37  is  approximately  the 
same  as  that  of  Figure  41,  as  expected. 


The  basic  conclusions  which  may  be  drawn  from  mode  formation  experiments 
in  this  C0?  laser  are  [49]: 

(1)  Far-field  burn  patterns  taken  at  1,000  and  520  nsec  are 
relatively  sharp  with  easily  seen  first  order  maxima.  Burns  taken  at  520, 

460,  440  nsec  show  a  successive  increase  in  loss  of  resolution  and  blurring. 

This  can  be  attributed  to  mode  formation  dynamics. 

(2)  This  pulsed  CO,  laser,  with  low  magnification,  unstable 

resonator  optics  in  a  decentered  configuration,  had  a  steady  state  transverse 
mode  formed  at  a  time  approximately  equal  to  that  predicted  by  the  Anan'ev 
formula  (Equation  1.2). 

(3)  A  steady-state  symmetric  transverse  mode  is  always  formed 
after  the  predicted  mode  formation  time.  No  competing  mode  of  similar  losses 
was  found  to  exist. 

D.  XeF  Laser  Experiments 

The  xenon  fluoride  laser  appeared  to  be  very  useful  for  inode  formation 
work.  The  initial  spontaneous  emission  was  much  higher  than  C0?  because  of  the 

8 

fourth  power  dependence  on  wavelength  -  a  factor  of  10  increase  in  initial 
intensity.  With  a  similar  average  gain  profile  the  XeF  laser  would  produce 
high  peak  power,  "turn  on”,  at  essentially  half  the  number  of  round  trips 
required  of  the  C02  laser  used  in  this  work.  Assuming  good  beam  quality,  the 

properties  of  the  mode  formation  should  therefore  be  readily  observed  early  in 
the  mode  formation  process.  This  ideal  approach  was  notably  modified  by 
experimental  problems  as  described  in  a  later  section. 

A  schematic  of  the  laser  used  is  shown  in  Figure  10.  The  laser  output 
aperture  is  10  x  20  cm  with  confocal  resonator  optics  of  separation  L  equal  to 
262  cm.  The  gain  length  is  200  cm.  No  discharge  voltage  is  used  in  the  XeF 
laser;  however,  a  pair  of  Helmholtz  coils  surround  the  cavity  to  supply  an 
e-beam  guide  field.  The  e-beam  voltage  is  considerably  higher  than  the  typical 
180  kV  gun  voltage  of  the  C02  laser,  normally  running  at  280  -  500  kV  and 

acting  as  the  sole  pump  mechanism  of  the  XeF  device.  The  e-beam  current 
through  the  foil  into  the  cavity  is  8  amps/cm2  current  density. 

1.  Gain  Measurements 

Before  installing  the  confocal  resonator  cavity  optics  on  the 
laser,  a  series  of  gain  measurements  was  made  by  determination  of  an  I/IG 

signal  with  a  short  XeF  pulse  laser  (20  nsec  FWHM)  as  a  function  of  high  energy 
pump  condition  (sidelight  fluorescence) .  The  sidelight  fluorescence  was  a 
measure  of  the  transverse  spontaneous  emission  of  XeF  at  0.3532  y,  thus  giving 
a  direct  reference  marker  of  small  signal  gain  in  the  laser.  Figure  44 
represents  a  moderate  number  of  shots  from  the  short  pulse  excimer  laser  for 
2  meter  (high  energy)  laser  pumping.  The  single  pass  I/IQ  signal  from  the 

small  laser  is  plotted  against  fluorescence  signal  from  the  2  meter  laser. 

This  figure  illustrates  the  point  that  threshold  gain  is  reached  essentially 
at  the  moment  of  sidelight  fluorescence  detection  for  the  M  *  1.32  magnification 
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FLUORESCENCE  F"  (mV) 


Figure  44.  Single  pass  I/Iq  versus  fluorescence  signal 

fluorescence  values  for  threshold  gain  are  shown 
on  the  curve  for  M  ■  1.32  and  M  *  2.7  resonators. 

resonator  optics  (0.8  mV)  and  a  short  time  after  this  value  for  the  M  ■  2.7 
resonator  optics  (3.6  mV).  The  purpose  of  the  gain  measurement  versus  side¬ 
light  fluorescence  is  to  show  that  the  onset  of  sidelight  fluorescence 
signifies  threshold  gain  in  this  laser,  within  10  nsec  or  less  for  the  M  *  1. 
resonator,  and  no  more  than  30  nsec  for  the  M  -  2.7  resonator.  This  may  be 
seen  from  Figure  45  which  is  an  oscilloscope  trace  of  sidelight  fluorescence 
with  the  resonator  mirrors  in  place  and  extraction  of  high  peak  power  (loaded 
cavity)  as  a  function  of  time.  To  determine  emission  of  high  peak  power 
relative  to  fluorescence  (threshold  gain),  the  temporal  profiles  of  these 
curves  were  electronically  summed  at  the  oscilloscope. 


Figure  45.  Sidelight  fluorescence  output. 

(Vertical  Scale  5  mV/div.,  horizontal 
scale  100  nsec/div.). 

2.  Far-Field  Measurements 

Figure  46  is  a  schematic  of  the  test  configuration  for  XeF  far- 
field  mode  formation  experiments.  A  helium  neon  alignment  laser  is  positioned 
such  that  the  beam  goes  through  a  2  mm  non-aluminized  centered  hole  in  the 
primary  mirror  and  strikes  the  secondary  (feedback)  mirror,  which  has  centered 
1  mm  reference  non-aluminized  spot.  Centering  the  helium-neon  beam  on  these 
holes  defines  the  optic  axis  of  the  resonator.  The  mirrors  are  aluminized  and 
the  He/Ne  alignment  beam  will  fill  the  laser  cavity  through  diffraction.  This 
resonator  is  then  aligned  by  centering  circular  interference  fringes  (formed 
between  successive  passes  of  the  laser  beam  in  the  resonator)  on  the  feedback 
mirror.  The  visible  beam  should  then  track  the  excimer  beam  and  can  be  used 
to  position  the  far-f ield  pattern  as  desired  (to  within  'v  50  microradians) . 

Because  of  external  optics  size  limitations  only  a  vignetted  portion  of 
the  laser  beam  was  used  in  far-field  mode  formation  analysis  for  both  the 
M  ■  1.32  and  M  ■  2.7  experiments.  In  the  case  of  the  lower  magnification 
resonator,  a  vertical  rectangle  of  1.21  cm  x  6.5  cm  was  extracted  from  the 
lower  side  of  the  beam  toward  the  anode.  For  the  M  5  2.7  resonator  a  small 
1.75  cm  x  0.34  cm  horizontal  rectangle,  centered  1.87  cm  below  the  feedback 
mirror,  was  the  extraction  aperture.  By  turning  flat  mirrors  these  beams  are 
propagated  approximately  29.6  meters  into  a  clean  room  area  to  an  off-axis 
parabolic  mirror,  10  cm  diameter  with  123  cm  focal  length.  All  mirrors  were 
characterized  to  X/10  or  better  at  0.3532  microns.  The  far-field  pattern  of 
the  off-axis  parabolic  mirror  was  then  magnified  by  a  5  cm  focal  length 
suprasil  lens,  which  later  imaged  the  pattern  at  over-exposed  developed 
Polaroid  film  planes  Cx  and  C2.  The  beam  was  first  attenuated  by  reflection 
from  a  fused  silica  wedge,  with  the  transmitted  portion  being  absorbed  on  an 
unexposed  developed  Polaroid  sheet.  The  reflected  part  of  the  beam  was  sent 
from  mirror  mx  to  a  30/30  chromium  beamsplitter  which,  after  suitable  reflec¬ 
tions  from  flat  mirrors  m2,  m3,  and  m4,  positioned  far-field  patterns  at  c2 
c2.  It  was  necessary  to  use  the  Polaroid  sheets  as  fluorescing  screens 
(exposed  Polaroid  type  52  sheets)  since  the  gated  pulse  cameras  were  insensitive 
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Figure  46.  Experimental  excimer  laser  configuration  for 
mode  formation  experiments. 

to  the  uv  output  from  the  laser  (0.3532  microns).  The  response  for  the 
Polaroid  sheets  was  remarkably  rapid.  The  sheets  were  checked  with  a  15  nsec 
pulse  from  the  gain  measurement  XeF  laser  and  no  pulse-broadening  from 
fluorescent  relaxation  could  be  detected. 

The  patterns  from  the  cards  at  Cj  and  c2  were  focused  on  two  high  speed 
image  converter  cameras.  Both  cameras  were  TRW  Model  10  image  converter 
cameras  equipped  with  three  exposure  framing  heads.  Each  of  these  cameras 
provide  three  pictures  with  arbitrary  exposure  times  of  5,  10,  or  20  nano¬ 
seconds  and  arbitrary  interframe  times  between  pictures  of  50,  100,  200,  and 
500  nsec.  The  five  nanosecond  exposure  time  was  used  in  all  of  the  frames 
taken  in  this  experiment.  This  was  done  to  minimize  the  far-field  pattern 
change  as  a  function  of  mode  formation.  One  of  the  fluorescing  cards  was  also 
viewed  by  a  S-5  vacuum  photodiode  (ITT  4018)  to  monitor  laser  pulse  length 


(not  shown  in  Figure  46).  Timing  reference  gates  from  the  framing  cameras 
were  electronically  added  to  the  photodiode  output  to  establish  frame  timing 
with  respect  to  the  laser  pulse.  At  the  laser  end,  vacuum  photodiodes  re¬ 
corded  laser  energy  as  well  as  cavity  sidelight  fluorescence. 

3.  The  Real  Versus  the  Ideal 

Three  serious  problems  were  encountered  during  these  tests: 
gas  life,  mirror  damage,  and  beam  steering.  The  last  two,  in  particular, 
tended  to  degrade  the  quality  of  the  experiments. 

Gas  life,  the  number  of  laser  pulses  with  a  useful  output/gas  fill,  was 
a  problem  with  the  stainless  steel  laser  chamber,  but  was  especially  bad  at 
the  beginning  of  these  experiments.  A  previous  experiment  with  this  laser  had 
utilized  XeCl  in  Raman  conversion.  Since  we  had  little  experience  with  the 
effects  of  residual  chlorine  on  the  XeF  mix,  the  mode  formation  study  was 
begun  with  only  a  perfunctory  cleaning  of  the  laser  chamber.  It  soon  became 
obvious  that  this  was  a  mistake,  as  we  were  getting  only  one  of  two  useful 
shots  per  laser  fill  as  compared  with  a  normal  six  to  ten  shots.  This  problem 
was  cured  by  a  total  disassembly  and  thorough  cleaning  of  all  laser  chamber 
parts . 

The  second,  and  more  serious,  problem  was  beam  steering.  The  steering 
requirements  were  far  higher  than  had  ever  been  demanded  of  this  system. 

The  initial  scheme  called  for  the  extraction  of  portion  of  the  laser  beam  n, 
1.21  x  6.5  cm,  steering  this  to  the  off-axis  parabola,  focusing  to  a  far-field 
pattern,  and  magnifying  this  pattern  so  that  the  central  five  diffraction 
angles  filled  the  10  cm  field  of  view  of  the  framing  cameras.  The  focal 
length  of  the  parabola  was  123.2  cm,  thus 


5fA  5(1.232)  (3.53.2  x  1(T9)  ,  , 

_  _____  !.8  x  10  m 

where  f  is  the  parabolic  focal  length  and  d,  the  slit  width. 
This  was  equated  to  the  camera  image  size 
1.8  x  lO-4  M  »  4  (0.0254) 

thus 

M  -  567. 


(2.9) 

(2.10) 

(2.11) 


With  this  system,  the  central  lobe  position  could  vary  by  about  - - 

and  still  be  in  the  field  of  view  of  the  camera.  This  required  the  laser 
beam  steering  to  be  repeatable  to  about  44  microradians.  The  initial  tests 
mentioned  above  showed  that  the  steering  did  not  meet  this  requirement.  A  set 
of  burns  was  taken  behind  the  lens  at  a  magnification  of  15  which  showed  the 
beam  on  a  shot-to-shot  basis  to  be  confined  to  an  area  about  two  Inches  in 
diameter.  This  implied  a  focal  point  movement  of  about  3  mm,  for  a  steering 
accuracy  of  only  2.4  milliradians.  This  gross  variation  was  not  understood 
since  the  laser  alignment  and  aiming  method  was  felt  to  be  far  better  than  the 
test  Indicated.  The  procedure  utilized  a  helium-neon  laser  on  the  laser 


optical  axis  which  was  transmitted  through  the  entire  system  to  aim  the  laser. 
The  resonator  was  aligned  by  interferometrically  overlaying  the  He-Ne  return 
spot  from  the  primary  and  secondary.  At  this  point,  the  laser  was  being  dis¬ 
assembled  for  cleaning,  so  it  was  decided  that  the  mirrors  should  be  aluminized 
to  allow  the  He-Ne  beam  to  fill  the  mode  volume  by  diffraction,  and  to  use 
the  resultant  Ne-Ne  beam  to  aim  the  laser.  Unfortunately  the  chlorine  in  the 
laser  chamber  had  already  reacted  with  the  dielectric  coatings  on  the  mirrors, 
and  when  cleaned,  large  portions  of  their  surfaces  showed  damage.  Due  to  a 
lack  of  time  (stripping  and  refiguring  would  take  several  weeks) ,  it  was 
decided  to  proceed  with  the  aluminizing.  The  M  *  2.7  optics  were  aluminized 
at  the  same  time.  The  resulting  mirrors  did  not  look  good,  but  appeared  to 
work  well  enough  for  the  experiment. 

With  the  aluminized  mirrors,  we  could  be  much  more  confident  in  our 
alignment  and  aiming  procedures.  The  aiming  problem  was  not  eliminated,  and 
it  was  determined  that  the  problem  was  due  to  a  combination  of  the  large 
secondary  mirror  mass  and  the  electron  beam  magnetic  guide  field  of  the 
laser.  When  the  laser  is  fired,  the  guide  field  coils  are  turned  on  about 
0.6  seconds  before  the  electron  gun.  The  mechanical  shock  to  the  system  was 
exciting  a  resonance  in  the  secondary  mirror  mounting  assembly  so  that  the 
mirror  vibrated  for  about  three  seconds  after  guide  field  initiation.  This 
was  observed  by  watching  the  He-Ne  spot  in  the  clean  room  and  firing  the  coils 
by  themselves.  The  problem  was  essentially  cured  by  leaving  the  coils  on  for 
3 %  seconds  before  firing  the  laser.  This  reduced  the  steering  error  to  that 
resulting  from  the  settling  of  the  secondary  mirror.  The  system  could  be 
aimed  with  an  accuracy  of  about  380  microradians.  In  later  experiments  with 
the  M  =  2.7  resonator,  the  mirror  resonance  with  the  magnetic  field  shocks  all 
but  disappeared,  and  an  aiming  accuracy  of  better  than  100  microradians  was 
achieved  with  the  former  magnetic  field  e-beam  delay  of  0.6  seconds.  The 
380  y-rad  steering  accuracy  was  not  sufficient  enough  for  the  experimental 
design  goal,  so  the  image  magnification  was  reduced  to  about  60  at  the  fluo¬ 
rescing  cards  to  ensure  obtaining  the  far-field  pattern  at  that  point. 

4.  Test  Results 

(a)  M  =  1.32  Resonator 

Because  of  the  degraded  mirror  quality,  the  far-field 
patterns  taken  by  the  image  converter  cameras  were  not  diffraction  limited. 

The  spot  sizes  were,  at  best,  about  five  times  the  diffraction  limited  central 
spot  size.  Even  with  this  restriction,  mode  formation  can  easily  be  seen  in 
Figure  47,  which  relates  mode  formation  time  to  XeF  far-field  burn  patterns. 
This  figure  shows  the  time  resolved  data  from  shot  number  3/5-2.  The  upper 
photo  is  the  fluorescence/laser  output  overlay.  The  center  photo  is  the 
photodiode  signal  in  the  test  room  overlaid  by  the  framing  camera  timing 
markers,  and  the  bottom  two  photos  are  the  framing  camera  pictures  of  the 
mode  evolution.  The  center  trace  determined  laser  threshold  relative  to 
emission  of  high  peak  power  seen  by  the  negative  photodiode  signal  with 
positive-going  framing  camera  gate  signals.  The  lower  photographs  were  taken 
at  image  planes  c1  and  c2  by  the  image  converter  cameras.  Frames  1  through  3 
show  images  which  were  "pinching  down,"  implying  the  mode  had  not  yet  formed. 
Frames  4  through  6  were  essentially  identical;  a  quasi-steady-state  mode 
had  formed.  From  these  data  it  is  apparent  that  the  mode  developed  between 
Frames  3  and  4,  for  a  mode  formation  time  of  between  320  and  430  nsec.  A 
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separate  run  showing  convergence  of  the  mode  along  the  lines  of  Figure  47  is 

shown  in  Figure  48.  In  this  case  mode  formation  occurred  between  360  and 

470  nsec.  In  both  these  figures  t  is  threshold,  t,  is  the  onset  of  high 

o  lase 

peak  power,  and  tcj£vert  c^e  crowbaring  of  the  e-gun,  ending  the  laser 

pulse.  The  far-field  pattern  in  Frames  4,  5,  and  6  is  a  blurred  area  with  a 
maximum  of  5X  diffraction  limited  dimension  and  is  representative  of  the  quasi 
steady-state  transverse  mode. 
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Figure  48.  Mode  formation  time  for  the  M  *  1.32  resonator 
(Shot  Number  3/5-1). 


Figure  49  is  a  summary  of  all  mode  formation  runs  taken  with  the  M  =  1.32 
resonator  from  the  vignetted  aperture  of  1.21  x  6.5  cm.  A  range  of  time  for 
mode  formation  is  shown  rather  than  a  specific  time.  This  is  the  result  of: 
(1)  poor  mirror  quality  and/or  the  XeF  active  medium  (since  the  far-field  spot 
size  was  at  least  five  times  the  central  lobe  size  of  a  diffraction  limited 
beam);  (2)  poor  resolution  of  the  framing  cameras  and  shadowing  of  images  by 
the  image  converter  tube  grid  wires  (these  cameras  represent  a  twenty  year 
old  technology);  and  (3)  differences  between  the  cameras  in  sensitivity  and 
actual  appearance.  Some  of  the  runs  in  Figure  49  are  open  ended — the  e-beam 
may  have  crowbarred  too  soon,  the  framing  intervals  may  have  been  too  widely 
spaced,  or  the  framing  intervals  may  have  been  incorrectly  timed,  giving  only 
a  lower  or  upper  bound  for  far-field  pattern  stabilization.  The  runs  were 
divided  almost  equally  between  a  mixture  of  1:3:760  NF3:Xe:Ar  at  two  atmos¬ 
pheres  and  1:3:1140  NF3:Xe:Ne  at  three  atmospheres.  Neon  diluent  appeared  to 
provide  a  better  beam  quality.  It  is  interesting  to  note  that  the  experi¬ 
mental  data  show  the  mode  to  have  formed  around  350  nsec — considerably  less 
than  the  Anan'ev  time  of  513  nsec.  This  anomaly  is  explained  in  detail  later 
in  this  work. 


(b)  M  =  2.7  Resonator 

Experiments  using  this  resonator  were  performed  in  a  manner 
similar  to  those  using  the  M  =  1.32  resonator.  The  resonator  optics  used  in 
this  test  series  were  in  better  condition  than  the  M  =  1.32  optics  and  gave  a 
somewhat  better  far-field  pattern.  The  short  mode  formation  time  for  this 
configuration  (168  nsec)  required  working  at  the  practical  limit  of  the  laser; 
i.e. ,  four  atmosphere  cavity  pressure  and  maximum  laser  gas  pumping.  The 
results  were  the  same  for  this  laser  as  the  upper  limit  of  the  C02  experiments. 

The  foil  separating  the  laser  and  the  e-gun  chambers  was  punctured  by  arcing 
after  a  few  runs. 


The  portion  of  the  beam  used  for  these  tests  was  a  0.34  cm  by  1.75  cm 
rectangle  located  below  the  feedback  mirror.  The  far-field  pattern  size  was 


increased  along  the  dimension  of  interest  by 


1.21 

.34 


3.6x  compared  to  the 


M  =  1.32  pattern.  Figure  46  shows  the  test  schematic.  To  compensate  for  the 
lower  transmitted  energy,  the  Suprasil  wedge  (B.S.l)  was  replaced  with  a  di¬ 
electric  coated  total  reflector.  A  typical  run  sequence  for  this  resonator 
is  shown  in  Figure  50,  which  is  a  counterpart  of  Figures  47  through  48  for 
the  M  =  1.32  case.  The  first  frame  is  not  resolved  because  of  the  low  in¬ 
tensity  from  the  laser  pulse.  It  appears  that  mode  formation  occurred  in  this 
run  between  130  and  170  nanoseconds. 


The  accumulated  data  for  the  M  =  2.7  resonator  are  shown  in  Figure  51. 

A  problem  of  a  dim  first  frame  applies  to  all  runs  from  shot  numbers  4/17-2 
through  4/17-5.  Again,  as  with  the  M  =  1.32  data,  these  data  show  an  anomalous 
early  mode  formation  of  time  of  140  nsec  as  compared  to  the  typical  theoretical 
result  of  168  nsec.  Data  were  taken  with  NF3:Xe:Ar  1:3:760  at  two  atmospheres 
and  NF3XE:Ne  1:3:1520  at  four  atmospheres. 


50 


fluorescence 

LASER  SIGNAL 
100  m.  /di*. 


LASER-MARKER 

SIGNAL 

100  m./div 


Figure  50.  Mode  formation  time  for  tiie  M  -  2.7  resonator 
(Shot  Number  4/17-3). 
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NOMINAL  MODE  FORMATION 
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Figure  51.  Mode  formation  in  the  M  =  2.7  resonator. 

A  summary  of  the  results  of  the  XeF  experiments  and  initial  modeling 
results  were  presented  at  the  Annual  Meeting  of  the  Optical  Society  of  America, 
at  Tucson,  Arizona  in  1982  [50].  These  results  will  be  submitted  for  publi¬ 
cation  at  a  later  date  [51]. 

5.  Concluding  Remarks  on  the  XeF  Laser  Experiments 

The  excimer  laser  was  anticipated  to  be  of  great  use  in  the 
verification  of  the  formulas  of  Zemskov  and  Anan'ev  for  transverse  mode 
formation  time  in  unstable  resonators  of  moderate  magnification  (1.32  <M  <  2.7) 
This  was  primarily  because  of  short  time  laser  "turn  on,"  i.e. ,  high  peak 
power  was  normally  achieved  well  within  the  Anan'ev  time  with  a  resonator 
magnification,  M  of  1.32  and  certainly  somewhat  within  the  time  for  the  higher 
magnification  of  2.7.  Therefore,  the  far-field  patterns  of  such  a  laser  could 
be  readily  examined  both  before  and  beyond  the  predicted  mode  formation  time. 

Various  experimental  limitations  did  not  permit  a  "fine  tuned" 
definition  of  mode  formation.  An  ideal  limit  would  have  produced  a  central 
lobe  of  the  far-field  pattern  rise  to  within  90  percent  of  its  peak  value,  and 
the  central  lobe  position  to  within  10  percent  of  the  nominal  angular  steady- 
state  value;  as  would  have  been  desired  to  validate  an  Anan'ev  criterion.  A 
much  coarser  stability  definition  came  from  these  experiments.  This  stability 
was  determined  as  variations  observed  from  a  nominal  blur  circle  five  times 
the  diffraction  limited  central  maximum  arising  from  a  severely  vignetted  laser 
aperture.  The  anomalous  shortening  of  mode  formation  time  will  be  justified 
in  Chapter  III  with  the  experimental  constraints  applied  at  this  juncture. 

Saturable  gain  was  suspected  to  contribute  to  a  large  degree  in  the 
observed  early  mode  formation.  The  excimer  laser  quickly  achieved  saturation 
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at  approximately  150  nsec.  Calculations  of  saturable  gain  for  this  laser 
(in  Section  III)  tend  to  support  the  hypothesis  of  early  mode  formation. 


The  early  mode  formation  time  has  been  substantially  confirmed  for  the 
XeF  laser  case,  taking  into  account  the  experimental  limitations  of  the  ex¬ 
periment  described  above  as  well  as  modeling  results.  A  series  of  experiments 
currently  being  formed  at  Maxwell  Laboratories  at  San  Diego,  where  these 
initial  results  were  taken,  may  confirm  this  general  assumption  [52], 

III.  THEORETICAL  MODELING 

A.  Introduction 

A  strip  resonator  called  CAV2D  (see  Appendixes)  has  been  used  to 
model  mode  formation  dynamics  [53].  This  has  been  developed  from  a  Fourier 
transform  code,  the  SOQ  code  of  Siegman  and  Sziklas  [54],  which  is  used  to 
determine  the  lowest  loss  steady  state  transverse  mode  in  unstable  resonators. 
The  code  can  model  external  vignetting  and  saturable  resonator  gain  using  the 
Rigrod  equations  [55, 56], in  an  infinitely  thin  slice  near  the  primary  mirror. 

A  fast  Fourier  transform  (FFT)  subroutine  [57]  has  been  incorporated  as  a 
considerable  time  saving  measure  for  the  large  Fresnel  number  cases  in  cal¬ 
culations  of  the  direct  and  inverse  Fourier  transforms.  As  an  example,  con¬ 
sider  a  typical  run  in  this  work  for  an  XeF  laser  with  a  tube  Fresnel  number 
:  2700.  16,384  mesh  points,  m,  were  used  in  the  calculation.  A  typical 

calculating  time  using  the  FFT  algorithm  is  proportional  to  m  log2  m  rather 
than  m2  for  the  usual  methods  of  propagation:  the  direct  solution  of  the 
wave  equation,  the  Huygens-Fresnel  integral  equation  [58],  or  the  Fourier 
transform  solution.  The  time  saved  is  enormous  for  large  Fresnel  number 
applications,  e.g. ,  the  ratio  of  calculating  time  is 
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(3.1) 


Considering  a  typical  time  of  300  seconds  for  one  using  the  FFT,  the  run 
by  alternative  computational  approaches  would  be  10  hours.  In  addition,  an 
expanding  coordinate  system  useful  for  the  unstable  resonator  geometry  is 
included  [ 59] . 

CAV2D  is  a  Fortran  program  which  calculates  the  transverse  modes  of 
unstable  resonators  for  lasers.  It  performs  wave-optics  calculations  for  the 
propagation  of  a  beam  back  and  forth  between  the  two  mirrors  of  the  resonator 
until  it  obtains  a  transverse  distribution  of  the  complex  optical  amplitude 
which  is  self-consistent,  1 . e . ,  which  repeats  itself  over  a  round  trip  to 
within  a  complex  multiplicative  constant.  The  outcoupling  of  a  strip  resona¬ 
tor  equals  unity  minus  the  square  of  the  absolute  value  of  the  multiplicative 
constant  (eigenvalue) .  This  is  the  obvious  purpose  of  the  code  for  steady- 
state  calculations  in  unstable  resonators. 

Typically,  a  starting  wave  which  has  uniform  intensity  and  has  a  phase 
distribution  equivalent  to  a  collimated  output  beam  is  launched.  This  wave  is 
assumed  to  have  just  been  reflected  from  the  convex  mirror.  The  starting 
wave  may  be  of  a  different  form  without  affecting  the  final  steady-state 
result.  One  possibility  would  be  a  wave  of  constant  intensity  with  randomized 
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phase  distribution,  though  this  choice  would  ordinarily  require  more  itera¬ 
tions  to  achieve  a  solution  than  would  the  uniform  phase  choice.  The  iterative 
method  can  be  justified  merely  as  a  numerical  means  of  obtaining  a  solution 
of  the  lowest-loss  mode,  without  any  reference  to  modeling  of  the  actual 
temporal  progression  of  the  course  of  events. 

Successive  propagation  back  and  forth  between  mirrors  is  the  physical 
nature  of  the  action  of  a  resonator.  Hence,  the  CAV2D  code  in  this  work  is 
extended  to  model  the  transverse  mode  formation  process  in  which  randomized 
radiation  from  spontaneous  emission  develops  into  a  well-formed  transverse 
mode.  The  code  extensions  involve  introducing  randomness  into  the  initial 
phase  distribution,  calculating  the  far-field  distribution  after  each  round 
trip,  saving  tables  of  values  of  1/1  and  the  directional  orientation  of  the 
point  of  maximum  far-field  intensity  after  each  round  trip,  and  printing  or 
punching  those  tables  for  further  statistical  analysis.  The  principal  features 
of  the  code  are  unchanged  from  those  ;ed  to  obtain  steady-state  mode 
distributions. 

The  code  specifically  treats  strip,  or  cylindrical-mirror,  resonators. 

The  limitation  of  one  transverse  dimension  is  not  as  restrictive  as  it  may  at 
first  seem.  Resonators  with  square  or  rectangular  cross-sections  have  the 
property  that  the  mode  distribution  for  an  empty  cavity  (equivalent  to  an 
unsaturated  uniform  gain  situation)  may  be  factored  into  a  product  of  a  strip 
resonator  mode  on  each  of  the  two  transverse  dimensions.  When  saturated  gain 
plays  a  significant  role,  or  when  the  symmetry  is  broken  (for  example  by  a  hole 
in  a  mirror),  the  factoring  is  no  longer  valid.  A  separate,  but  somewhat 
similar,  computer  program  called  CAV3D  has  been  prepared  and  used  for  the 
cases  in  which  there  is  a  hole  in  the  concave  mirror.  This  code,  which  takes 
into  account  two  transverse  dimensions,  is  capable  of  handling  Fresnel  number 
cases  in  the  realm  of  CC^  resonators  (N^,  <  100)  but  is  certainly  not  valid  in 

the  regime  of  excimer  lasers  (2700  <  <  11,000)  because  of  computer  memory 

limitations.  Even  one  of  the  largest  computers,  e . g. ,  the  CDC  7600  with  the 
inclusion  of  auxiliary  "large  core"  memory  (LCM) ,  does  not  have  the  capability 
of  handling  256  x  256  mesh  points. 

A  verification  of  CAV2D  for  steady-state  transverse  modes  has  been 
achieved  by  comparison  with  the  results  of  other  workers  using  an  entirely 
different  methodology.  In  particular,  the  asympotic  expansion  code  results  of 
Horwitz  [16,  21]  and  the  work  of  Chen  and  Felsen  [60]  were  examined.  The 
quite  appreciable  detail  which  is  present  in  near-field  intensity  distributions, 
readily  apparent  from  the  inspection  of  plots  in  the  Horwitz  papers,  are 
accurately  reproduced  by  the  CAV2D  code. 

The  angular  distribution  of  the  far-field  intensity  (1(6)]  is  calculated 
in  CAV2D  after  each  round  trip  in  the  resonator  by  direct  numerical  evaluation 
of  the  Fraunhofer  integral  written  as 
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where  0  is  a  diffraction  angle  and  C  is  a  normalization  constant. 


Considering  one  printout  per  round  trip,  the  bulk  of  paper  generated  over 
one  mode  formation  run  requires  condensation.  The  criterion  for  data  selection 
becomes  two  independent  data  points  per  round  trip  selected  as  follows: 

(1)  the  intensity  of  the  -highest  peak  in  the  far-field  distribution 
is  normalized  to  an  intensity  for  a  beam  with  identical  power,  uniform  near¬ 
field  intensity  and  phase  across  an  aperture  of  width  D; 

(2)  the  angular  position  of  the  highest  peak  is  normalized  to  the 
diffraction  angle  (X/D)  of  the  aperture. 

The  total  number  of  data  points  depends  upon  the  number  of  round  trips 
acceptable  for  convergence  of  the  data  to  a  steady-state  condition.  Thirty 
to  forty  round  trips  per  run  appear  to  be  suitable  for  C02  and  XeF  with 

M  =  1.32.  No  more  than  twenty  round  trips  per  run  are  acceptable  convergence 
criterion  for  XeF  with  a  resonator  magnification,  M,  of  2.7. 

B.  C02  Laser  Analytical  Modeling 

A  saturable  gain  model  for  the  C02  experiments  has  been  built  into  the 
code  based  upon  C02  test  parameters.  A  nominal  average  gain  (0.009/cm)  and  an 

assumed  initial  intensity  at  threshold  intensity  relative  to  saturated  in¬ 
tensity,  I/I  of  5  x  10-12  was  used.  "hese  initial  conditions  were  deter¬ 
mined  by  experimental  measurement  of  average  small  signal  gain  during  the 
period  of  mode  formation;  and  measurement  of  loaded  cavity  gain,  y,  relative 
to  small  signal  gain,  y  ,  at  21  round  trips  after  threshold  in  the  laser  at 

the  approximate  mode  formation  time.  The  ratio  I/I  can  be  derived  at  this 

sac 

time  from  the  ratio  of  loaded  cavity  gain  to  small  signal  gain 
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Yq  is  obtained  from  Figure  20  from  the  small  signal  intensity  at  approximately 
21  round  trips  (I 'T  7.5)  and  from  gain  length,  5,; 
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y  is  derived  from  the  average  increase  in  intensity  of  one  round  trip  during 
the  gain  switched  spike  leading  edge  at  approximately  21  round  trips 
(I'/l  '  ^  3.5),  as  seen  in  Figure  23,  so  that 
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Therefore,  using  Equations  (3.4),  (3.5),  and  (3.6),  we  have 
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which  determines  the  degree  of  saturation  in  the  actual  lasing  process  for  a 
specific  number  of  round  trips.  This  information  specifies  the  initial 
conditions  of  the  analytical  code  saturable  gain  inputs  for  the  C02  laser 
cases. 

Figure  52  shows  one  selected  run  for  the  C02  case  with  saturable  gain 

included.  One  value  of  relative  normalized  intensity,  I/Iq,  and  angular 

position  of  maximum  peak  in  terms  of  relative  diffraction  angle  D0/X  is  shown 
for  each  round  trip.  Thirty  round  trips  are  shown  in  the  figure.  The  nominal 
convergence  of  the  graphs  is  around  the  mode  formation  time  as  predicted  by 
Anan'ev  for  22  round  trips.  The  data  presented  in  Figure  52,  when  combined 
with  19  other  runs,  give  cumulative  statistical  results  on  mode  formation  in 
the  C0„  laser  with  saturable  gain  conditions  as  shown  on  later  figures  as 

normalized  to  a  mode  formation  time  of  20  round  trips. 

The  convergence  of  the  far-field  distribution  as  a  function  of  round 
trip  number  is  typically  shown  by  Figures  53  through  62.  Again,  each  round 
trip  peak  intensity  has  been  normalized  to  unity  from  a  laser  aperture  of 

dimensions,  D,  plotted  on  a  scale  of  ±  —  where  A  is  the  laser  wavelength. 

A  random  phase  wavefront  was  launched  with  initial  conditions  similar  to 
those  of  Figure  52:  I/I  at  threshold  is  5  x  10“12  and  average  small 

Sat 

signal  gain,  A  ,  is  0.9%/cm.  The  variation  of  the  far-field  shape  and  position 

is  shown  in  the  series  of  figures  to  converge  to  the  steady-state  far-field 
shape  and  position  is  shown  in  the  series  of  figures  to  converge  to  the 
steady-state  and  far-field  pattern  (Figure  62) . 

The  compilation  of  20  runs  with  differing  values  of  random  phase  but 

with  the  same  initial  conditions  as  that  of  Figure  52  is  shown  in  Figures  63 

and  64  for  the  relative  intensity  of  the  highest  peak,  I/I  ,  and  the  relative 

0D  ^ 

diffraction  angle,  —  ,  of  the  highest  peak,  respectively.  These  are  shown 

statistically  with  the  mean  (M)  and  one  standard  deviation  (S)  for  each  round 
trip.  These  tend  to  converge  to  the  typical  relative  Anan'ev  number 
(Equation  1.1)  of  20  round  trips  (22  round  trips  in  reality).  A  line  is  drawn 
across  the  statistical  compilations  which  represents  the  mode  formation  round 
trip  number  as  determined  from  the  Anan'ev  formula  (Equation  1.1).  This 
represents  a  baseline  for  mode  formation  and  is  referred  to  as  the  "nominal 
mode  formation"  in  all  figures  of  this  type.  The  effect  of  gain  introduced 
in  these  calculations  seems  to  produce  a  change  in  apparent  mode  formation 
time.  This  may  be  seen  from  a  single  run  comparison,  for  example,  of  Figure  65 
and  Figure  52.  In  Figure  65,  the  wavefront  is  launched  with  the  same  initial 
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Figure  56.  Far-field  pattern  at  4  round  trips. 
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Figure  59.  Far-field  pattern  at  9  round  trips. 
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Figure  62.  Far-field  pattern  at  24  round  trips. 

random  phase  numbers  as  in  Figure  52  but  without  gain.  A  comparison  of 
Figures  52  and  65  for  the  first  few  round  trips  is  as  expected — small  gain  in 
the  loaded  cavity  gives  results  comparable  to  that  of  the  empty  cavity  since 
intensity  is  not  depleting  the  gain  (small  signal  gain  region).  However, 
with  the  onset  of  saturated  gain  at  18  round  trips  the  situation  begins  to 
change.  The  relative  intensity  of  the  dip  located  at  approximately  21  round 
trips  is  much  reduced  in  Figure  52  than  in  Figure  65.  This  causes  a  con¬ 
vergence  at  22  round  trips  in  Figure  52  instead  of  27  as  in  Figure  65.  The 
angular  convergence  in  both  figures  is  relatively  unaffected.  The  overall 
effect  appears  to  be  a  somewhat  earlier  mode  formation  time. 

The  statistical  accumulation  of  the  20  empty  cavity  runs  with  the  identi¬ 
cal  initial  phase  conditions  of  runs  generated  for  Figures  63  and  63  is  shown 
by  Figures  66  and  67,  respectively.  The  initial  conclusion  of  earlier  mode 
formation  time  in  the  loaded  versus  unloaded  cavity  cases  is  borne  out  by 
the  statistical  analysis  of  Figures  66  and  67. 

The  particularly  poor  convergence  of  the  l/lD  statistics  of  the  empty 

cavity  runs  as  compared  to  the  loaded  cavity  cases  is  apparent  in  Figures  63 
and  66.  An  inspection  of  the  individual  runs  of  Figure  66  reveals  that  certain 
runs  did  not  converge  as  expected  but  continued  to  change  well  beyond  the 
Anan'ev  criterion.  This  appeared  to  be  caused  by  oscillation  in  the  relative 
intensity  values,  and  could  be  interpreted  as  competition  of  two  or  more  modes 
of  equal  loss  which  would  by  necessity  lengthen  mode  formation  time.  The 
application  of  typical  gain  in  the  C02  case  seems  to  minimize  if  not  eliminate 
this  degeneracy.  This  potentially  explains  why  mode  competition  has  not  been 
observed  at  mode  crossing  points  which  are  located  at 

F  -  n  +  0.875  .  (3-8) 

eq 


63 


RELATIVE  ROUND  TRIPS 


RELATIVE  ROUND  TRIPS 


279. 


INTENSITY  (l/l0) 


NUMBER  OF  ROUND  TRIFS 


Figure  65.  Empty  cavity  calculation  of  mode  formation  under 
identical  starting  phase  conditions  of  Figure  52 
without  gain,  M  =  1.32,  N  *  279. 


!  Figure  66.  I/Iq  statistics  of  20  runs  for  CO^  with  empty  cavity, 

M  -  1.32,  N  -  279. 
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We  chose  a  Fresnel  number  for  these  calculations  to  produce  maximum  mode 
separation.  However,  the  low  magnification  minimizes  mode  separation  and 
emphasizes  mode  competition.  This  could  be  inherent  within  the  definition  of 
"mode  formation  time"  for  increasingly  smaller  magnification,  especially  in 
laser  systems  which  do  not  saturate.  In  short,  saturable  gain  appears  to 
clear  up  any  potential  mode  competition.  Angular  deviations  in  Figures  64  and 
67  did  not  particularly  appear  to  be  affected  by  gain. 

In  one  series  of  experiments,  mode  competition  in  the  CO^  laser  was 

deliberately  sought.  This  was  done  by  spacing  the  resonator  mirrors  to  give 
the  "bad"  equivalent  Fresnel  number  of  Equation  (3.8).  The  two  competing 
modes  were  calculated  in  empty  cavity  codes  to  have  basically  the  same  far- 
field  pattern  but  differing  phase  angles.  This  should  have  given  a  temporal 
beat  frequency  easily  observed  in  the  near-field  laser  pulse  with  a  photon 
drag  detector.  However,  no  such  beat  frequency  was  observed  and  far-field 
burn  patterns  were  not  abnormal.  The  tentative  conclusion  is  that  gain  must 
eliminate  or  at  least  minimize  the  effects  of  mode  competition  in  loaded 
resonators  compared  to  the  empty  cavity  results  calculated  by  other  authors 
at  mode  crossing  points  [61].  The  formation,  of  a  transverse  mode  at  an 
earlier  time  with  saturable  gain  for  the  empty  cavity  modeling  results,  may 
be  analyzed  in  a  similar  manner  as  stable  resonator  modes  [62].  In  the  empty 
cavity  results,  the  lowest  loss  mode  is  obtained  through  diffraction  losses  in 
the  resonator.  With  saturable  gain,  a  nonlinear  effect  is  added.  The  lowest 
loss  mode  tends  to  predominate  earlier — by  reducing  the  gain  available  to 
other  competing  modes — with  the  result  that  the  mode  forms  earlier  than  with 
the  empty  cavity  condition. 

It  will  be  shown  for  later  cases  of  the  excimer  laser  modeling  at 
M  =  1.32  that  saturation  seen  early  in  the  pulse  at  10  round  trips  (compared  to 
30  round  trips  for  mode  formation)  readily  gives  an  earlier  mode  formation  for 
loaded  cavity  than  for  empty  cavity  calculations  and  is  well  within  the 
Anan'ev  time. 

C.  XeF  Laser  Analytical  Modeling 

This  work  was  similar  to  the  CO  modeling  but  with  the  requirement  of 

the  larger  tube  Fresnel  numbers  of  2700  for  the  M  -  1.32  resonator  magnifica¬ 
tion  case  and  10,800  for  the  M  =  2.7  case.  a  addition,  the  initial  spontaneous 
intensity  was  much  higher  (a,  1  x  10-4)  for  I/l  -  in  the  excimer  case  at 

threshold).  The  average  gain  was  fortuitously  siir.iaar  to  that  of  C02  at  0.9%/ 

cm.  Saturation  was  achieved  under  these  conditions  much  earlier  than  in  the 
CO  cases  at  approximately  1/3  of  the  XeF  mode  formation  time  for  M  ■  1.32  at 

10  round  trips. 

Another  restriction  was  the  aperture  vignetting  necessary  in  the  excimer 
laser  cases.  The  M  ■  1.32  runs  were  analyzed  through  a  near-field  vertical 
rectangular  stop  of  1.21  cm  width  by  a  height  of  6.5  cm.  The  2.7  runs  were 
even  more  restricted  by  a  horizontal  rectangle  of  0.34  cm  width  by  1.75  cm 
length.  These  tended  to  change  the  mode  formation  time,  apparently  by  a 
spatial  filtering  process. 
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In  these  cases,  the  tube  Fresnel  number  was  considerably  larger  than  in 
the  C02  runs  (279  versus  2,700  -  10,800).  The  C02  runs  were  valid  at  1024 

mesh  points:  16,384  were  adjudged  necessary  for  the  XeF  cases.  The  10,800 
tube  Fresnel  number  was  not  as  critical,  in  any  event,  since  experimental 
results  were  not  as  precise  for  the  rapid  convergence  of  data  in  that  magnifi¬ 
cation  regime. 

The  statistical  compilation  of  15  runs  with  the  gain  conditions  of  the 
M  -  1.32  XeF  laser  are  shown  in  Figures  68  and  69.  These  cases  include 
external  vignetting  with  a  1.21  x  6.5  cm  rectangle  in  the  laser  beam.  Mode 
formation  has  been  obviously  moved  to  shorter  times  with  respect  to  the 
Anan'ev  formula.  This  appears  to  be  partly  a  result  of  the  vignetting  process 
but  primarily  of  early  gain  saturation.  Saturation  was  achieved  in  Figures  68 
and  69  after  approximately  7  relative  round  trips  (10  actual  round  trips) . 
Relative  round  trips  are  defined  in  the  statistical  runs  as  the  Anan'ev 
time  normalized  to  20  round  trips. 

A  comparison  must  be  made  of  these  data  with  those  of  Figures  70  and  71 
for  empty  cavity  calculations  with  identical  initial  random  phase  starting 
conditions  for  each  run.  Again,  as  in  the  C02  cases,  saturable  gain  plays  a 

part  in  mode  formation.  In  the  loaded  cavity  (Figure  68)  the  mode  is  formed 
after  approximately  14  relative  round  trips,  whereas  in  the  empty  cavity 
(Figure  70)  it  has  developed  after  approximately  18  relative  round  trips. 

The  angular  development  of  the  mode  is  not  appreciably  different  in 
either  case  (Figures  69  and  71) .  Vignetting  appears  to  move  mode  formation 
to  slightly  earlier  times  compared  to  the  relative  Anan'ev  time. 

The  empty  cavity  results  have  been  statistically  analyzed  for  the 
M  =  2.7  cm  case  in  Figures  72  and  73.  Thirteen  runs  similar  to  those  of 
Figures  70  and  71  were  utilized  in  the  calculations.  These  were  run  at  a 
mesh  point  dimension  of  16,384.  One  individual  run  was  obtained  as  a  compari¬ 
son  to  Figures  72  through  73  with  32,768  mesh  points.  The  convergence  results 
were  similar.  This  represents  the  maximum  memory  capability  of  the  CDC  7600 
computer  for  these  calculations.  An  accurate  reproduction  of  the  near-field 
phase  and  intensity  for  the  M  =  2.7  resonator  would  require  approximately 
60,000  mesh  points.  The  consequence  of  a  smaller  number  of  mesh  points  is  a 
smoothing  of  the  near-field  phase  and  intensity,  giving  a  spatially  filtered 
effect  in  the  far-field.  This  appears  as  a  secondary  influence  on  mode 
formation. 

Vignetting  is  more  important  in  the  M  *  2.7  case  because  of  the  smaller 
vignetted  aperture  which  causes  more  spatial  filtering.  Excluding  any  gain 
effects,  the  mode  formation  time  has  been  moved  back  to  16  round  trips  as 
opposed  to  20. 

Finally,  in  this  series,  consider  Figures  74  through  77  which  were  run 
with  full  aperture  for  the  excimer  case  of  M  ■  1.32,  ■  2700.  These  cases 

were  run  for  comparison  with  previous  vignetted  examples  to  show  earlier  mode 
formation  from  vignetting  and  also  to  directly  evaluate  saturable  gain  effects 
in  large  aperture  excimer  lasers.  The  mode  formation  time  was  found  to  move 
forward  from  14  round  trips  (Figure  68)  to  16  round  trips  (Figure  74)  presumably 
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Figure  72.  I/Io  satistics  of  13  runs  for  XeF  empty 
calculations  and  vignetting,  M  -  2.7, 
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because  of  the  absence  of  vignetting  effects.  Also,  the  angular  selection  is 
not  as  rapid  for  the  nonvignetted  case  (Figure  75)  as  for  the  vignetted  case 
(Figure  69).  As  a  concluding  example.  Figures  It  through  77  show  full  aper¬ 
ture  from  runs  for  the  excimer  laser  case  M  =  1.32,  NT  =  2700,  without  gain. 

A  comparison  of  these  figures  to  Figures  74  through  75  shows  that  saturable 
gain  does  significantly  affect  mode  formation  time.  The  same  initial  random 
phase  of  the  runs  in  Figures  68  through  69  was  used  to  generate  Figures  74 
through  75  and  Figures  76  through  77  for  direct  comparison. 

In  general,  early  mode  formation  time  of  the  XeF  experimental  data  may 
be  explained  by  Figure  68.  Without  respect  to  experimental  considerations 
of  poor  far-field  resolution,  the  mode  has  formed  early  in  time.  In  fact, 
the  figure  shows  mode  formation  at  about  14  relative  round  trips  (about  360  nsec) 
rather  than  20  relative  round  trips  (513  nsec) .  An  examination  of  the  ex¬ 
perimental  results  for  the  M  =  1.32  resonator  of  Figure  49  gives  a  comparable 
answer.  Therefore,  the  modeling  and  the  experiment  agree  in  terms  of  the 
results.  Mostly  as  a  result  of  saturable  gain,  however,  the  mode  is  formed 
ahead  of  the  nominal  Anan'ev  time. 

A  consideration  of  the  M  =  2.7  results  produces  a  similar  conclusion. 

The  extreme  vignetting  in  this  case  moves  the  mode  formation  to  an  earlier 
time  than  it  is  in  the  M  =  1.32  case.  The  empty  cavity  calculations  of 
Figures  72  through  73  show  the  mode  to  have  been  formed  around  16  relative 
round  trips.  This  is  a  time  of  134  nsec  which  is  consistent  with  the  results 
of  Figure  51.  In  this  case,  because  of  larger  magnification,  mode  competition 
is  not  a  serious  problem,  and  with  earlier  mode  formation  saturable  gain  should 
not  be  as  much  of  an  effect  in  moving  back  mode  formation  time.  Spatial 
filtering,  because  of  extreme  vignetting,  should  be  the  predominating  variable. 

IV.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  geometrical  formulation,  as  proposed  by  Anan'ev  for  the  formation  of 
transverse  modes  in  positive  branch  unstable  resonators  as  an  upper  limit,  has 
been  validated  experimentally  and  by  statistical  modeling  for  moderate 
resonator  magnifications  (1.32  £  M  2.7)  over  a  large  regime  of  tube  Fresnel 
numbers  (279  <  N^  <  10,800) . 

A.  Summary  of  Experimental  Results 

The  objective  of  these  experiments  was  to  observe  the  far-field 
pattern  from  two  selected  high  energy  lasers,  CO2  amd  XeF.  at  times  before, 

during,  and  beyond  nominal  transverse  mode  formation  time,  as  predicted  from 
geometric  optics  theory.  This  represents  the  transition  time  required  for  the 
laser  to  convert  the  radiation  in  the  cavity  from  initial  spontaneous  emission 
to  a  quasi-steady-state  self-replicating  wavefront  characteristic  of  a  dif¬ 
fraction  limited  laser  beam. 

The  initial  temporal  points  of  the  lasers  were  referenced  to  threshold 
gain  conditions  experimentally  determined  by  small  signal  gain  measurements. 

Mode  formation  was  determined  from  the  threshold  time  to  the  time  required 
for  steady-state  formation  of  the  far-field  pattern. 


78 


RELATIVE  ROUND  TRIPS 


The  CO^  experiments  examined  far-field  burn  patterns  of  truncated  pulses 

which  were  varied  before  and  beyond  the  nominal  mode  formation  time  by  varia¬ 
tion  of  laser  gain.  These  data  indicated  mode  formation  in  the  range  of 
460-520  nsec,  agreeing  well  with  the  predicted  value  of  490  nsec  for  an 
M  »  1.32  resonator. 

The  XeF  experiments  used  an  electronically  imaged-gated  framing  camera 
to  observe  a  fluorescing  image  in  the  far-field.  The  advantage  of  the  XeF 
laser  was  that  it  "turns  on"  early  with  high  peak  power,  allowing  the  far- 
field  to  be  observed  in  one  series  of  experiments  at  approximately  half  the 
predicted  mode  formation  time.  These  data  indicated  mode  formation  in  the 
regime  of  300  -  400  nsec  (predicted  mode  formation  time  513  nsec)  for  an 
M  =  1.32  resonator  and  130  -  180  nsec  (predicted  mode  formation  time  168  nsec) 
for  an  M  =  2.7  resonator. 

The  apparently  anomalous  early  mode  formation  time  of  the  XeF  cases  can 
be  justified  by  modeling  analysis  of  saturable  gain  and  laser  aperture 
vignetting,  as  well  as  by  poor  far-field  patterns.  Poor  far-field  patterns 
were  the  result  of  severely  degraded  resonator  optics  and  poor  imaging  quality 
of  framing  cameras.  These  characteristics  made  it  difficult  to  determine  a 
reasonably  exact  time  of  formation  of  the  mode. 

The  results  of  the  C02  work  indicate  a  formation  time  consistent  with  the 
Anan'ev  prediction.  The  XeF  work  gives  a  significantly  earlier  time  for  mode 
formation. 


B.  Summary  of  Theoretical  Modeling 

The  CAV2D  code  used  in  this  work  was  a  modification  of  earlier  codes 
which  derived  the  lowest  loss  transverse  mode  in  unstable  resonators.  It 
used  Fourier  propagation  within  the  resonator  to  develop  a  self-consistent 
steady-state  wavefront  output. 

The  code  had  some  time  saving  devices  built  into  it:  the  fast  Fourier 
transform  (FFT)  algorithm  and  an  expanding  coordinate  system. 

Instead  of  being  interested  in  the  steady-state  transverse  mode,  we 
desired  to  obtain  wavefront  properties  of  phase  and  amplitude  at  each  round 
trip  in  the  resonator  at  an  early  time  and  project  these  properties  to  the 
far-field  for  statistical  analysis.  A  random  phase,  uniform  amplitude  wave 
was  launched  to  begin  a  convergence  process  as  a  function  of  resonator  round 
trip . 


Statistical  modeling  of  the  C02  case  showed  convergence  around  the 

Anan'ev  time  and  agreed  with  the  CO^  experiments.  Modeling  with  vignetting 

in  the  XeF  cases  showed  the  convergence  time  of  the  lowest  order  steady- 
state  mode  to  be  considerably  earlier  than  the  Anan'ev  mode  formation  time. 
This  was  consistent  with  the  experimental  evidence. 

An  interesting  result  of  the  analytical  work  was  that  saturable  gain 
early  in  time  had  an  appreciable  effect  in  developing  a  steady-state  pattern 
in  the  far-field  considerably  sooner  than  in  the  empty  cavity  case  or  Anan'ev 
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predictions.  With  the  knowledge  of  early  saturable  gain  in  certain  high 
energy  lasers,  it  appears  that  the  Anan'ev  formula  is  quite  conservative  and 
that  the  actual  transverse  mode  forms  earlier  in  time. 


C.  Recommendations  for  Future  Work 

It  would  be  desirable  to  repeat  the  XeF  experiments  under  better 
conditions  which  would  allow  more  accurately  determined  mode  formation  times 
in  the  range  of  resonator  magnifications  (1.3  <  M  <  3)  of  interest.  The 
modifications  suggested  for  a  more  accurate  set  of  results  are: 

(1)  the  use  of  better  resonator  optics  (better  than  A/20  RMS 
surface  figure) , 

(2)  the  use  of  technologically  current  gated-image  cameras 
having  more  frames  and  greater  spatial  resolution  per  run),  and 

(3)  the  use  of  the  full  aperture  of  the  laser  for  far-field 
analysis,  which  would  necessarily  require  large,  good  quality  optics  to  direct 
a  10  x  20  cm  beam  (the  size  of  the  "2  Meter  Maxwell  laser"  beam) . 

Recent  communications  with  individuals  at  Maxwell  Laboratories, 
Incorporated,  at  San  Diego,  California  indicate  that  work  in  these  directions 
is  now  in  progress  [52]. 
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CAV2D  CODE  DESCRIPTION  [63] 


CAV2D  is  primarily  a  two-curved-mirror  code,  which  is  realistic  for  most 
unstable  resonators.  The  most  common  mirror  configuration  for  unstable 
resonators  includes  a  concave  mirror  which  is  at  least  as  large  in  transverse 
extent  as  the  gain  region,  and  a  convex  mirror  which  is  smaller  than  the  gain 
region.  A  substantial  fraction  (the  outcoupling  fraction)  of  the  radiation  is 
propagated  from  the  concave  mirror  toward  the  convex  mirror,  passing  outside 
of  the  convex  mirror  and  is  emitted  fror  the  resonator.  The  convex  mirror  is 
taken  as  mirror  number  one  in  the  code,  while  the  concave  mirror  is  mirror 
number  two.  For  convenience  the  calculation  is  normally  started  with  a  wave 
which  has  just  been  reflected  from  the  convex  mirror.  For  ordinary  calcula¬ 
tions  not  involving  modeling  of  mode  formation,  this  starting  wave  is  uniform 
in  intensity.  It  also  has  a  phase  distribution  which  is  the  same  as  would 
apply  to  an  idealized,  i.e. ,  geometric-optics,  resonator  wave. 

The  most  common  configuration  is  a  confocal  or  near-confocal  one.  However, 
this  limitation  is  not  imposed  in  the  code.  Instead,  the  code  accepts  as 
input  data  the  g  values  (XG1  and  XG2)  of  each  mirror.  Near  statement  10  there 
is  a  calculation  of  the  normalized  (i.e . ,  measured  in  units  of  the  mirror 
separation)  radius  of  curvature  (RP1  and  RP2)  of  each  of  the  mirrors,  and  the 
corresponding  magnifications  (XMAG1  and  XMAG2)  associated  with  propagation 
from  mirror  1  to  mirror  2  and  from  mirror  2  to  mirror  1,  respectively.  When 
the  input  g  values  correspond  to  a  confocal  resonator,  the  magnification 
XMAG2  is  unity,  corresponding  to  a  collimated  beam  propagating  from  mirror  2 
toward  mirror  1,  while  the  magnification  XMAG1  equals  the  overall  magnifica¬ 
tion  which  in  magnitude  equals  the  ratio  of  radii  of  curvature  of  the  two 
mirrors.  The  advantage  of  dealing  with  g  values  is  that  these  quite  general 
formulas  are  based  on  g  values  as  independent  variables.  Of  course,  the  code 
could  start  with  input  radii  of  curvatures  and  calculate  the  g  values. 

Near  statement  8  the  code  tests  the  g  values  to  see  whether  they  corres¬ 
pond  to  the  unstable  resonator.  If  the  resonator  is  stable,  the  code  rejects 
that  set  of  input  data,  prints  a  message  (format  86)  "STABLE  RESONATOR,"  and 
proceeds  to  the  next  set  of  input  data.  If  desired,  this  could  be  modified  to 
permit  calculations  of  stable  resonators.  The  code  detects  whether  the  mirrors 
are  reversed  from  the  normal  ordering  in  which  mirror  1  is  the  convex  mirror. 

If  the  ordering  is  reversed,  a  message  is  printed  (format  96)  "BACKWARD  WAVE 
CALCULATION"  and  the  values  of  XMAG1  and  XMAG2  are  suitably  adjusted.  This 
feature  has  only  been  used  for  limited  tests. 

A  factor,  named  FACT,  which  is  the  ratio  of  the  number  of  total  mesh 
points  (XNTOT  or  NTOT)  to  the  number  of  mesh  points  across  each  mirror  (XMIR 
or  XNMIR) ,  enters  into  various  steps.  The  total  Fresnel  numbers  XNF1T  and 
XNF2T,  which  are  involved  in  the  propagation  kernals,  are  greater  than  the 
input  values,  XNF1  and  XN2T,  of  the  ordinary  Fresnel  numbers  of  the  mirrors 
by  the  square  of  this  factor. 

The  ratio  of  the  transverse  size  of  the  mirrors,  A10A2,  and  its 
reciprocal,  A20A1,  are  needed  in  determining  the  location  of  points  involved 
in  the  interpolation  process.  They  are  obtained  as  the  square  root  of  the  ratio 
of  Fresnel  numbers  of  the  two  mirrors,  and  its  reciprocal.  These  are  cal¬ 
culated  below  statement  3. 


If  the  two  g  values  are  equal  and  the  two  Fresnel  numbers  are  equal,  the 
resonator  is  symmetric  (this  special  case  is  often  used  as  discussed  else¬ 
where)  .  [23]  A  test  for  this  is  made  at  consecutive  statements  126  and  112. 
If  the  resonator  is  symmetric,  a  message  (format  83)  "SYMMETRIC  RESONATOR"  is 
printed,  and  a  simple  formula  is  used  to  calculate  the  equivalent  Fresnel 
number  FEQ.  It  it  is  not  symmetric,  a  message  (format  84)  "NONSYMMETRIC 
RESONATOR"  is  printed;  and  FEQ  is  calculated  by  a  more  involved  formula.  The 
latter  is  based  on  the  formulation  for  an  equivalent  symmetric  resonator 
(which  assumes  that  the  actual  concave  mirror  is  oversized  as  compared  to  the 
beam) .  GEQ  is  the  g  parameter  of  the  "symmetric-equivalent"  resonator,  and 
XNSEQ  is  its  Fresnel  number.  The  relevant  formulas  are  given  in  the  "crib 
sheet."  In  either  event,  the  value  of  FEQ  is  printed.  The  special  feature, 
which  treats  a  ring  resonator  with  unequal  lengths  of  the  two  paths,  will  not 
be  considered  here. 

To  determine  when  the  iterations  have  converged,  use  is  made  of  the 
ratio  of  the  complex  amplitude  at  a  test  position  on  two  successive  round 
trips.  This  position  is  determined  by  the  value  of  an  input  parameter  TSTFR 
(whose  default  value  is  zero) ,  such  that  the  position  is  removed  from  beam 
center  by  an  amount  equal  to  TSTFR  times  the  half-width  of  the  overall  beam 
(including  guard  band) . 

CAV2D  Program  Modules 

The  program  consists  of  the  following  modules,  of  which  all  but  the  MAIN 
program  are  subroutines: 


MAIN 

GGNPM 

PROPCBM 

GGUBS 

FOURT 

BQUAL 

INTRP 

PHSEM 

NFLPP 

FRFLD1 

GAINM2 

MFLPP 

The  four  modules  listed  in  the  right-hand  column  are  of  relatively  minor 
signf icance  but  have  been  listed  for  completeness.  The  GGNPM  and  GGUBS  sub¬ 
routines  are  taken  from  the  IMSL  Library  which  is  available  on  most  mainframe 
computers.  If  RNPH2  is  non-zero,  a  correlated  randomness  of  phase  at  the 
position  of  mirror  2  is  applied  to  simulate  effects  such  as  gas  density 
fluctuations.  The  random  numbers  for  that  purpose  are  obtained  from  GGNPM, 
but  this  is  not  of  much  current  interest.  For  mode-formation  calculations,  a 
randomness  is  applied  to  the  phase  of  the  starting  wave.  This  is  called,  if 
the  parameter  RANST  (mnemonic  for  "random  start")  is  non-zero.  The  pseudo¬ 
random  numbers  used  for  that  purpose  are  obtained  from  the  subroutine  GGUBS. 
BQUAL  calculates  a  simple  measure  of  beam  quality,  but  it  is  not  as  useful  as 
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calculating  the  overall  far-field  intensity  distribution.  PHSEM  can  be  used 
to  introduce  a  "phase  modulation"  across  the  beam,  but  is  probably  not  needed 
for  future  work. 

The  MAIN  program  controls  the  overall  calculations.  It  reads  input  data, 
initializes  various  arrays,  registers  and  prepares  the  "starting  wave;" 
then  begins  the  iterative  process  of  calculating  propagation  back  and  forth 
between  the  two  mirrors.  This  iterative  process  is  continued  until  it  either 
satisfies  a  convergence  criterion  (based  on  input  or  default  value  of  CVTST 
parameter)  or  completes  the  maximum  prescribed  number  of  iterations,  XMITR. 

The  program  always  prints  the  values  of  data  parameters,  which  determine  the 
problem  being  calculated,  and  generates  at  least  one  line  of  print  for  each 
iteration,  which  gives  an  indication  of  the  iterative  behavior.  It  optionally 
calls  the  line  printer  plotting  program  NFLPP,  and  produces  plots  of  intensity 
and/or  phase  across  the  entire  aperture.  An  earlier  version  could  optionally 
plot  across  one-half  of  the  mirror;  but  the  current  version  does  not  permit 
this,  as  it  is  not  appropriate  if  the  resonator  is  not  ideally  aligned. 

Control  parameters  PLLP1  and  PLLP2  control  the  plots  at  position  1  and  2. 

Such  plots  may  optionally  be  made  only  after  completion  of  the  iterations  (if 
PLLP1  or  PLLP2  equals  unity)  or  after  each  iteration  (if  PLLP1  or  PLLP2  equals 
two) .  Upon  completion  of  one  problem,  the  MAIN  program  reads  input  data  for 
another  problem,  if  such  is  available. 

The  overall  program  is  arranged  such  that  the  typical  phase  distribution 
is  an  idealized  flat  distribution  across  the  reference  surface,  i.e. ,  a 
pseudo -collimated  beam  is  being  propagated.  Of  course,  there  are  departures 
from  such  ideal  behavior  at  the  starting  point  for  each  propagation  step,  and 
departures  that  arise  from  diffraction  as  a  result  of  each  propagation  step. 

The  MAIN  program  establishes  the  complex  amplitude  distribution  U(I)  prior  to 
each  propagation  step,  then  calls  the  PRPCBM  (mnemonic  for  "propagate  colli¬ 
mated  beam")  subroutine.  The  MAIN  program  also  furnishes  to  the  PRPCBM  sub¬ 
routine  the  value  of  the  overall  Fresnel  number  (called  XNFT  in  the  PRPCBM 
subroutine)  for  the  propagation  step,  and  a  value  of  the  parameter  XPAVG  as 
well  as  the  total  number  of  mesh  points,  NTOT.  XPAVG  is  most  commonly  set 
to  unity.  If  XPAVG  is  greater  than  unity,  the  propagation  kernel  is  averaged 
over  XPAVG  subintervals  at  each  mesh  point.  The  PRPCBM  subroutine,  in  turn, 
calls  the  fast-Fourier-transform  subroutine  FOURT.  The  values  of  complex 
amplitude  after  propagation  are  returned  in  the  same  array  U(I).  FOURT  is 
the  only  subroutine  that  is  called  by  another  subroutine;  all  others  are  called 
by  MAIN. 

After  each  return  from  the  propagation  subroutine  PRPCBM,  the  MAIN  pro¬ 
gram  can  optionally  call  the  plotter  program  NFLPP.  In  addition,  after  re¬ 
turn  from  PRPCBM  for  propagating  from  mirror  2  back  to  mirror  1  (i.e.,  back  to 
the  position  of  the  convex  mirror) ,  there  is  an  optional  call  to  FRFLD1  for 
calculation  of  the  far-field  intensity  distribution.  The  call  to  FRFLD1  is 
made  for  each  iteration,  if  the  control  parameter  CFFLD  is  greater  than 
unity;  this  option  is  useful  in  mode  formation  calculations.  Aside  from 
these  optional  calculations,  the  interpolation  subroutine  INTRP  is  called 
after  completion  of  a  propagation  step  in  order  to  obtain  the  complex  ampli¬ 
tude  distribution  across  the  "receiving"  mirror  in  a  suitable  form.  The 
propagation  steps  have  generated  the  complex  amplitude  distribution  at  points 
which  cannot  cover  the  receiving  mirror  with  NMIR  mesh  points,  as  desired.  For 
example,  in  the  propagation  from  concave  to  convex  mirror  of  a  confocal 
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resonator,  the  calculated  distribution  ha.',  the  range  of  NMIR  points,  covering 
the  full  size  of  the  beam;  whereas  one  wants  to  determine  a  new  distribution 
with  NMIR  points,  covering  the  small  size  of  the  convex  mirror. 

After  interpolating  to  find  the  distribution  across  the  "receiving" 
mirror,  there  is  an  optional  call  to  the  plotter  program  NFLPP  to  plot  the 
distribution  across  the  entire  aperture,  as  referred  to  earlier. 

After  completion  of  the  two  propagation  steps  and  associated  options,  the 
program  saves  the  ratio  of  the  amplitudes  (at  a  selected  position)  for  the 
last  two  iterations.  This  ratio  is  used  to  generate  a  test  parameter  TEST, 
which  approaches  zero  as  the  calculation  converges,  i.e.,  as  the  ratios  for 
successive  round  trips  tends  to  become  the  same  with  each  new  round  trip. 

The  absolute  value  of  the  difference  from  unity,  called  ABTS,  is  tested  again 
in  a  preselected  convergence  test  parameter,  CVTST.  When  it  is  less  than 
CVTST  on  two  successive  evaluations,  the  calculation  is  assumed  to  be 
successfully  converged  (as  signalled  by  the  parameter  STOPTS  being  equal  to 
two) . 


After  convergence,  there  is  an  optional  plot  of  near-field  distribution 
as  noted  above.  Also,  there  is  an  optional  calculation  of  the  far-field 
intensity  distribution  by  the  FRFLD1  subroutine.  Often,  the  far-field  in¬ 
tensity  distribution  is  also  plotted.  For  mode  formation  calculations,  there 
would  be  excessive  printed  output,  if  the  far-field  distribution  were  plotted 
each  time  it  is  calculated,  namely,  after  each  iteration.  Hence,  an  option 
is  provided  to  omit  the  far-field  pattern.  For  mode  formation  studies,  the 
values  of  IMXFF(ITER)  and  of  FFMAX(ITER)  are  calculated  in  the  FRFLD1  sub¬ 
routine  and  returned  to  the  MAIN  program.  For  simplicity,  these  values  are 
always  generated,  even  if  CFFLD  is  not  greater  than  unity,  so  as  to  produce 
mode  formation  plots  as  discussed  below.  These  are  the  values  of  the  angular 
position  of  the  most  intense  far-field  point  and  the  relative  intensity 
(compared  to  the  value  for  a  uniphase  near -field  output  with  the  same  in¬ 
tegrated  intensity)  at  that  point,  for  the  present  iteration,  as  indicated  by 
the  index  ITER. 

Upon  completion  of  the  iterative  process,  the  MAIN  program  may  optionally 
(if  CFFLD  is  greater  than  unity)  call  the  "mode  formation"  subroutine  MFLPP. 

The  subroutine  produces  a  line-printer  plot  of  the  values  of  FFMAX(ITER)  and 
IMXFF(ITER)  for  inspection.  The  arrays  are  printed  with  identifying  parameters 
to  logical  unit  number  1  for  punching  into  cards  for  further  analysis.  This 
writing  to  unit  1  will  be  suppressed  if  PCHMF  is  less  than  or  equal  to  zero. 

After  completion  of  one  problem,  control  is  transferred  to  the  beginning 
of  the  program  (by  a  GO  TO  1  statement)  to  read  input  data  for  a  new  problem. 
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APPENDIX  B 
CAV2D  CODE  LISTING 


CAV2D  CODE  LISTING 


DIMENSION  U!Jt7t>8),  UP ! 32  766 ) 

DIMENSION  RNORM!32768),CRAND!I276e) 

dimension  i*xrr (2jo ,rrMAx(iOu> 

COMPLEX  U,UP 

COMPLEX  RAT  1 .SAvUl .SAVU1P 

COMMON  /COM1/u/CON2/UP/COn3/RNO«m/CGM4/CRAND 
LEVEL  2,  U  .UP  »  RNORM  ,  C  R  A  ND 
COMPLEX  PRMT.CFKTfi 

c. this  oouble-precison  needed  et  nature  of  subroutine 

DOUBLE  PRECISION  ISEED.JSEED 
61  FOPMAT!1X,4F10.5.c10.4.F12.5) 

7C  F0RMAT!1HC.9HC0NVERGED> 

71  FGRMAT!1Hu.20HMAX  NO  OF  ITERATIONS) 

76  FORMATHHu) 

77  FORMAT(IHI) 

7 •  FORMATCiX.lMHERROR  IN  CAVT7) 

83  format!'  symmetric  resonator') 

64  FORMAT!'  NON-SYMHETRIC  RESONATOR') 

£5  FORMAT!'  EQUIVALENT  FRESNEL  NUMBER  FEO=',F1C.5> 

56  FORMAT!1X,1tHSTABLE  RESONATOR) 

37  FORMAT!'  SPECIAL  TREATMENT  OF  RING  RESONATOR  WITH  UNEQUAL  PATHS') 
88  FORMAT!'  OF  COLLIMATED  ANO  EXPANDING  BEAM'S.  F£OP  =  ',F1C.S> 

39  FCRMATI1Xf4MX&l*,Fl0.5,5X,4HXG2*,Flu.5,5x,5MXNF1«,F10.5, 
15x,5mxnF2«,F1u.5,5X,5hEP$L*,F10.$> 

9?  FORMAT!1X,6HXNT0T*,f1C.5,5X,6HXNMIR«,F1'j.5) 

9 1  FORMAT!1X,6HXMAG1*,F10.5,5X,6MXMAG2*,F1*.5) 

92  FORMAT  !1X,6hXNF1T*,F1L.S.5X,6hXNF:t*,fK.5) 

93  FORMAT!'  XMITR  RAT1AM  KATIPH  &AM.RT  AVFLX  OUT 

TCP') 

94  FORMAT!'  EMPTY  RESONATOR  CALCULATION') 

95  FORMAT!'  LOADED  RESONATOR  CALCULATION  WITH  GAIN  IN  NARROW  REGION 
1 NE A h  'tHRCR  2.  GAIN  SATURATION  EY  RIGROD-LIAE  SCHEME.') 

NAMELljT/CVDAT/XGl.XGt.XNFI.XNFt.XNTOT.XNITR.XLAST. 

1XNTOR.CFFLD,  XNMIR.PA1I2,  XPAVG,  EPSL.GL,  RANST, 

i  TVBQ.aDOLP.  PLLP1.PLLP2. 

Z  PH AMI  .PHAM2  ,PHCTl  ,  PHC Y 2  ,GNO F S  ,  1NCPL  .BMTLT  ,  X I  NT S • 

*.RNPH2 ,ISEED,NCRLT,JSEED,VRCRN,XSPAM,XMM2F  ,  X F C TR , AMP R F  ,PHS R F , I S* M F , 


C**»*****TSTFR  ADDED  26  NOV  60. 
C  **•••**  *R  TPMS  ADDED  27  JAN  61. 


5 STOPN.XLPOL.TSTFP.PRDGS.CVTST .xinint.fr  inf .UMTlN.RTPHS.PCHMF, 
6FNFCS,SPCLFF,FRCMSK 
DATA  STOPN/2.C7 
DATA  XG1.XG2/1.25.1.25/ 

DATA  XNF1.XNF2/12.8.12.8/ 

OATA  XNT0T.XNMIR/512.0.256.C/ 

OATA  XN1TR/1  .0/ .XLAST/C.  :/ .XNTCR/3.U/ ,  CFFLD70.07 

DATA  GL/O.O/,  PAl  12 / 2 • j/ i  EPSL/0.0/.XPAVG/1.CF 

DATA  TVBQ/:.:/,uDOLP/0.-/,rMTLT/C.C/ 

DATA  PLLP1/..U/.PLLP2/C.O/ 

OATA  TSTFR/C.O/ 

OATA  PHAM1/ J.<7fPHAM2/C.0/,PMCYT/»,.0/,PNCr2/?.C/,GNOFS/0.O/ 

DATA  INCPL/1/  .X1NTS/1  .C/,RANST/u.0F 

DATA  RNPMt/0.0/,JSEED/55555.C/,t.CRLT/12/,ISlED/J3JJ1.0/,NRCRN/i:.0/ 
DATA  XSPAM/.,*/ 

DATA  XHM2F/«.L/,XFCTR/C.C/,AMPRF/u.0/,PHSRF/C.C/ 
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DATA  JSXPf/1/,PCHMF/C.C/ 

DATA  XLP0l/1.l/,PR0GS/C.C/ 

OATA  PI/J. 141592654/ 

DATA  CVTST/j.wCI/ 

DATA  XIMNT/O.OtO  /  ,  f  A  I  N  f  /  0.  C  .  Cl  1  /  ,fc  MT  I  N  /  .  .  „  / 

DATA  kTPHS/».*/,FNfCS/C.O/,SPCLFF/O.0/,fRCMSN/.''.n/ 

NX*  3276? 

C***»»REAO  AND  wRITE  INPUT  OATA  AND  DO  INITIAL  C AL C UL AT  1 0 NS 
1  READ  <5,CV0AT) 

If  <XLAST)9?,4,9l01 
A  wRIT£(C,CVOAT) 

STP  *  STOfN-1.0 
NTOT  *XNTOT 

9999  IMNTOT  .LE.  NX)GO  TO  9  5  C 
WRITE (6 ,9AC) 

940  f ORPAT  <  *  INVALID  RUN  SINCE  NTOT  WAS  GREATER  THAN  NX.') 

STOP 

950  CONTINUE 

If  (XNITR  .LE.  2C0.OC0  TO  952 
WRITE<6,951) 

951  FORMAT!'  XNITR  EXCEEDS  2Cj.0  NOT  SUITABLE  Rt  IMXFF  AND  FFMAX. 
STOP 

952  CONTINUE 
NNIR*XNMR 

RAT1TS*1.C 
STOPT.C.O 
XPITR*1 .0 
G1G2*XGl»xG2 

C...«..*STABLE  RESONATORS  are  rejected. 

If  (G1G2)K,«.6 

8  If !G1G2-1 .09,9,12 

9  WRITE(6,?6> 

GC  TO  1 

10  XAPG*G1G2*!g1G2-1.C> 

XARG*SORT ! X  A  RG ) 

I F ( G 1 G  2  .LT.  S.l)XARG*-XARG 
RPl*X62/ <XASG*G1G2“XCt) 

RP2*XGl/(XARG«G1G2-XG1) 

XMAG1*1.C*1.C/RP1 
xPAG2»1.0*1«0/RP2 
If !xnf2-xnF1)  3,5,5 
3  wRITE(6,9t) 

96  FORMAT!'  BACKWARD  WAVE  CALCULATION') 

XSAVV*XMAG1 
XKAG1*1 .C/XNAG2 
XKAG2*1 .C/XSAVV 
5  FACT-XNTOT/XNMIR 
FACT«f ACT**< 

XNF1T*XNF1»FACT*XMA01 

XNF2T»XNF2*FACT*XMAG2 

xnf;t*xnf2T/xlpol 

AAG*XNf i/XNf 1 
A20A1«SWRT(ARG) 

AlOA2*1.:/AiOAl 
WRITE (6,76) 

WRITE (6,89)XG1  ,XG.,XNF1  ,XNf2,tPSL 
WAITE !6,9l>XNT0T,XNPIR 


! 
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1  AGE 
•  F2T 


>2  > 
r*6^) 

••2-1 .G) 

3)/(xrAG1*XLP0L> 
WRITE  <6,87) 
.RITE<6,88)FEUP 


2.C*<TSTFR*1 .0) 


GO  TC  12^ 

CTR 

•  XhP2F>/2.r-  *1.3 


SRF  ) 

RNT) 

IRROR  has  central  APERTURE 


L  X  HUU 

t,'  LIHUU*',I6) 


»ITm  COMPLEX  REFLECT1 


OES. •••*•• 

/2 . 0*1 •  u 

.3 

TOT.RNORH) 


1 


pMarwin  j 

I* 


DO  15  J»HPLJ,LI*UJ 
15  $aand*$rand*pnorn<j ) 

C  R  A  N  D ( 2 ) “SRAND/SQRT (  X  F 

IMwacrn  .he.  C.O)  «r; 
20  CONTINUE 
122  DO  13b  1*1  ,nTOT 

U  <I)*CI1PLX<0. 0,0.0) 
UP<I>*Cnplx(0. 0,0.0) 
138  CONTINUE 

ARG*SQRT(XInTS> 

DO  113  I*LIPL.LINU 
U  <I)*C«PLX(AAG,0.0) 
110  CONTINUE 

rSEEO*ISEEO 
IFCRANSTJIO?.,  105,102 
132  DO  10A  1*LIML«L1MU 
J*I-LINL 

I F ( J  .LE.  10  )WRIT£(6, 
137  foRPAT(1x,2e30.20) 

CALL  GGUBSUSEED.I  .RAN 
ARG*(2.3*RAND-1.0)»PI* 
U  (I)*U  (1)*CEXP(CPLX 

10A  continue 

IF(XSPAP)11S,103,118 
118  NNID*XP1D 

DO  117  I*l1NL,NMID 
LLL*NT0T-I*1 
U  (I)  *  ( U  (I)  ♦  XSPAP 

U  ClLL)  *  XSPAN  •  U  (I 
117  CONTINUE 

1 C  3  IF(dPTLT)136,109,106 
106  DO  105  I*LIHL,LI*U 
XI*I 

ARG*(Xl-XMD)/XN«IR*2. 
U  <I)*CEXP<C*PLX<C.3,A 
105  CONTINUE 

109  IF(XG1-XGE)121  ,120,121 

i;o  ofsti*epsl*xnnir/2.o*( 

0FST2*0FST1 
GO  TO  128 

121  xnAGRT*XPAG1*a*AG2 

OF  ST1 *£PSl*xnhI R/2 .0*1 
OFST2*0.0 

C. ............  START  A  round 

1 2 S  CONTINUE 
SLl A**0.C 

DO  129  1*LIPL,LI»U 
SulA“*SulA»*UlI)*CCNJL 
129  CONTINUE 

SAVUl*U(ITEST) 

IF(Pha«1  ,NE.  0.3JCALL 
CALL  PROPCB (U ,NT0T , XNF 
IFCCPLLP2  .gt.  *.E>  .a 
1  iSTOPT  .GE.  STP)))  ca 
c*.. ••INTERPOLATE  TO  GET  U2P 
CALL  INTRPCU.U  P 

:••••••••« ig*od-ttpe  gain  at 
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AD-A136  515 
UNCLASSIFIED 


TRANSVERSE  MODE  FORMATION  IN  POSITIVE  BRANCH  UNSTABLE 
RESONATORSIU)  ARMY  MISSILE  COMMAND  REDSTONE  ARSENAL  AL 
DIRECTED  ENERGY  DIRE..  R  W  JONES  AUG  83 

DRSMI/RH-83-4-TR  SBI-AD-E950  467  F/G  20/5  NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  Of  STANDARDS-1963-A 


i 


IF(GL>131.1JS,131 

131  CALL  GAINR2CU  P ,NTOT ,N*IR ,Gl .GDGLP .GNOFS) 

£•*•«. PR0PA6ATE  to  HIMOI  1. *•••.*•**•••••**•...*•*•**••**•*••*. ..*••••• 
135  lF<RNPH2)33t36f33 
31  00  35  ISL1NL.L1MU 
AR6aCRAN0(X)*RNPM2 
U  P(1)«U  P(1)*CEXP<CHPLX(3.0.ARG)) 

35  CONTINUE 

36  CONTINUE 
SU2PAN*0»Q 

DO  141  I>LlMLiLlMU 

141  SU2PAN*$u2PAN*UP«X)*CONJGIUPEX)) 

X0EN*LIRU>L1HL*1 

avplx*$u2pam/xden 

00  143  1*1 »NTOT 

IFCPROGS  .GT.  C.3)VRITE<6,2Ei1)I,UP(I) 

2a1  FORMATE'  AT  140  UPE 16 , ') »' ,2E2: . 12 ) 

140  UEI)»U  P(I) 

1MPHAK2  .RE.  0.3KALL  PMSEP  EPMAMi  ,PHCT2  ,NPXR  ,NTOT  ,U  ) 

XINAHP*$QRT EXIN1NT) 

1FEXHH2F  .EQ.  0.9)  GO  TO  1044 

IF ELIMLL  .GT.  LIMUU)  GO  TO  1044 

DO  1043  I*LINLL.LINUU 

UE1)-UE1)*CFKTR 

IFEXININT  .EQ.  0.0)60  TO  1343 

XI«I 

ARC><XI-XKI0)/XNH1M2.C 

OLPH  *PI*FR1NF*ARG**2  ♦  EXMITR-1.0)  *  RTPHS 
OLPHI  a  OLPH  ♦  CXX-XHXO)/XHHXR*i.j«PI«BHTIN 
AR6C*C0$E0LPH1) 

ARGS*3 IN  EDLPH1 ) 

CARP»CHPLX<ARGC,-ARGS) 

UINCR*XINAMP*CARP 
WMTE(6,13tO)l»UINCR,DLPH,DLPHI 
1003  FORHATE1X,I4,4E20.12) 

U(I)  *  U ( X )  ♦  UXNCR 

1043  CONTINUE 

1044  CONTINUE 

1945  CALL  PR0PCBEU,NT0T,XNF2T,XPAV6) 

1M((PLLP1.GT.  9.3)  .AND.  ((XNXTR  .EQ.  XPITR)  .OR. 

1  (STOPT  .GE.  STP  )  )  )  .0R.EPLLP1  .GT . 1 .«) )CALL  NFtPP E J  ,  1 ,NTOT ,P A ll 2, 
21 . 1NCPL. XNXTR) 

IF (CFFLD  .GT.  1.0)CALL  FRFLOKU  ,NTOT ,NHIR  , EPSl , XMA6 2 . A1 0A2 , 
IXMXTR.INXFF.FFMAX.CFFLO.FNFCS'SPClFF.FRCMSK) 

C*#***INT£RPOLATE  TO  GET  U  P.  I.E.  AMPLITUDE  RECEIVE#  »T  MIAROR  1, 

CALL  INTRPEU.U  P ,LI«L ,LIPU ,OFST 1 , A10A2 .XNTOR ,XHAG 2) 

£.*... .GET  POaER  RATIOS***** 

Su1PAH*C.3 

00  145  1*L1RL«LIHU 

145  Sill  P4Ma$Ul  PAN*U  PEI)*CONJGEU  PCD) 

(UTAH  a  C.9 
OC  144  I»1 .NTOT 

144  SUTAH  a  (UTAH  ♦  UEX)  •  CCNJGEGCD) 

OUTCP  a  1,?  -ESUlPAM/SUTAM  •  AlO«2> 

GAMRTaSUlPAH/SUI AH 
SAVUlPauP(ITEST) 

C*****6£T  COMPLEX  RATIOS  OF  AMPLlTUOES  AT  MIRROR  CENTERS. ••••••*••*••••• 
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R«T1*$AVUlP/SAVUl 

UTIMUnoaii) 

UTItMCMSdATI) 

WRITE(6,61 )AN2TR,AATlAN,RAT1PM,GANRT,AVf LA.OUTCP 

TESTaRATlAH/AATITS-1.0 

RATITSMATlAN 

tlTtMBS(TUT) 

IF<A6TJ-CVT$T)159,159t1$8 

158  STOPT>0.0 
60  10  16C 

159  STOPTaSTOPT‘1.0 

160  CONTINUE 

. . ‘CALCULATED  VALUE  09  TV86L  XS  RATIO  01  CENTRAL  FAR-flELR 

{••••••INTENSITY  ARXS1N6  f RON  GEONETRIC'OPTICS  PORTION  Of  6EAN  TO 

{••••••VALUE  aMXCH  WOULD  RESULT  If  OUTPUT  MERE  UNIPHASE < 

If ( (TVBQ  .GT.  0.0)  .AND.  ((ANITA  .EQ.  ANITA)  .OR. 

KSTOPT  .GT.  1.0)))  CALL  BQUAL(NTOT,NNIR,U) 

IMSTOPT-  STOPN >169,190.190 

169  If (RN1TRoRN:TA)99,  191,170 

170  KNITR>ANITR«1.0 

C*»***RENEW  U  (I)  FOR  NEAT  ITERATION. ••«•••••*•*••*•*•**•••••• 

1M6D174, 173,17* 

173  BO  172  I-1.NT0T 

U  (I)»U  P(I)/RAfl  . 

If (PROGS  .GT.  0.0)URITE(6,202) I ,UPfI> 

202  fORNAT('  AT  172  UP(',I6,')-',2E20.1t) 

172  CONTINUE 
GO  TO  176 

17*  00  175  1-1.NT0T 
U  (I)-U  P(l> 

175  CONTINUE 

176  GO  TO  128 

190  WRITE (6,70) 

GO  TO  192 

191  WRITE (6,71) 

192  CONTINUE 

OPTIONAL  fAR-f IELO  CALCULATION ••*•••••*••••••••**•*••*• •*•••*••••• 

If  (CFfLO  .GT.  O.OCALL  FRfLBl  (U  ,NTOT  ,NNIR  ,EPSl,APA62,A10A2, 

1 ANITA, lNAff,fFNAA,CffLB,fNPCS,SPCLFf,FRCNSK) 

If  (CFfLO  *6T.  1.0)  CALL  Nf  LPPdNAPF  ,  f  FNA  A,  ANITA,  ISANf  , 
1ANF2,APA61,EPSL,ANNIR,TSEED,PCNNF) 

GO  TO  1 

99  WRITE(6,78) 

9001  STOP 
END 

SUBROUTINE  PR0PCB(U,NTCT,ANfT,APAV6) 

{••••DOES  COLLINATED  BEAN  PROPAGATION  Of  COPPlEA  ARRAT  U,  NXTN  NTOT 
(••••NESH  POINTS.  ANf T  IS  TOTAL  fRESNEL  NUNBER .  APAVG  IS  NUNBIR  Of 
{••••POINTS  USED  IN  AVERAGING  Of  PROPAGATION  KCANfcL (NONINALLV  UNITV). 
LEVEL  2,  U 

DIMENSION  WOAfc(IOO)  ,NN(1 > 

OINENSXON  U(NTOT) 

COPPLEA  U.ffTC 
BATA  PI/ 3. 1*1 59265 36/ 

NN(1)>NT0T 

RNTOT«NTOT 


XhTOT-NTOT n 
NPAVC-XPAVC 

FACTO. C/(l6.WI*XNFT) 

CALL  F0UAT<U.NN,1.«1,1,U0RX) 

00  132  X-1.NT0T 
«X»I 

fftc»chplx<:.o,c.o> 

00  120  J-1.NPAVS 
K2JM1>2*J-1 

XIX*XHT0T-XX«1.5-X2JH1/2.0/XPAV(. 

XASe»(XHTOT-ABStXX !>)••£ 

PSO*4.0*PI**2*XASO 

U(MSI*FACI 

ARCC»C0E<AR6> 

AI(S*tIH(AI() 

12C  FFTC-FFTC*CMPlX<AR6C.-AR6S> 

mcMFic/mvt 

132  U<I >*U<I >»F FTC 

CALL  F0URT<U.FiN.1.-1.1«W0RK> 

00  134  1*1  .NTOT 
134  U<1)«U<1)/XNT0T 
RETURN 
CNR 

SUBROUTINE  INTRP  CU .UP  ,LINL ,LIHU  ,OFST  , AROAS.XNTOR.XPAC  » 
C***BESSEL  FORMULA  FOR  INTERPOLATION  BV  HALVES.  SCARBORO  P  S3  ER  VI 
LEVEL  2.  U.UP  ( 

COMPLEX  U<1).UP(1>.CXNT1»CXNT2,CINT3.VAL 
ARC*ABS(XMA6> 

RTPAS-SSRTXARC) 

XHXb>(FL0AT(LlHL)*FL0AT(LlMU))/2.C 

00  130  I«LZML*LIMU 

ll*t 

ARCOXI-XNXD-OFST  i/xmac  *aroas«xmio 

1ARS-ARC 

XIARC-IARS 

FR ACT ■ARC- XI ARC 

VFRAC-FRACT-9.S 

CINT1*U<ZARC«1)-U(XARC> 

VAL*(U(IARG)«U(1ARC*1>)/2.0«C1NT1*VFRAC 

IF (XNTOR-2. 0)139.130.137 

137  CINT3«U(1ARS*2)-3.0*U<IARC«1)*3.0*U(IARS>-U(IARC-1) 
VAL*VAL«CXM3/6.0*VFRAC«(VFRAC**2-Q.2»> 

ISO  ClNT2*<U<IARC«2)-U<XARCO>-UtlARC>»U(lAAC-1>>/2.C 
VAL«VAL«CINT2/2.9*(VFRAC**2-C.2S> 

139  UP  (I)>VAL 
130  UP  (X)«UP  (D/RTHAC 
RETURN 
IkO 

SUBROUTINE  PNSEMXPNAM.PmCT 
LEVEL  2.  U 

COMPLEX  UCnTOD.CVAR  .PMSE 
OATR  PI/3 .1413920330/ 

XNT0T*NT0T 
XNMXR>NMIR 

XMXO>CXfcTOT«1.?)/2.3 
00  117  XO.NTOT 
XI-I 


> 


XAM6  •(XI«XRID>/XNH1R*PI*PMCV  *2.L 
TARO  *PH  AH  * (COS (XARG  )-1.0) 

CTAA  ■CHPLX(C.O,VARG  I 
PHSE  «C£XP(CTAR  > 

U(!)«U(1)*PHS£ 

117  CONTINUE 
RETURN 
end 

SUBROUTINE  NFLPP(CU,NMIN,NMAX,PA1I2, IPOS, 1NCR, XMITR > 


LEVEL  2,  CU 

DIMENSION  6RAPH (1C1 ) 
COMPLEX  CCENH.CARG 

COMPLEX  CL(NMAX) 

76  FORMAT (1H0) 

77  FORMAT  M Hi ) 

79  FORM A T (64 X, IMPOSITION, I 3) 
60  FORMAT (19Xi3H0.C,32X, 'XXX 

RELATIVE 

U.5 

INTENSITY 

XXX', 

131X.3H1.C) 

81  FORHAT(19x,3HO.0,52x,'xxx 

RELATIVE 

i.S 

AMPLITUDE 

XXX', 

131X, 3h1 . 0) 

85  FORMAT(18X,SH-PI/2,33X,'«« 

RELATIVE 

0  6  0 

PHASE 

36X,SH*PI/2> 

89  F0RHAT(63X,9MITERATI0N,I4) 
111  FORMATUOX.IOlAl) 


112  FORMAT(13X,l5,2X,13lAl) 

DATA  PI/3. 141562634/ 

ORTA  BLANK/IN  / , DOT/ 1 H ./, X/ 1MX/ , ZERO/1 MO/ . PLU S/1H»/ 
N6RAPh*101 
XNGRM1*nSRAPH-1 
••••••PRINT  HEADING 

WRITE (6 .77) 

ITCR-XPITR 
WRITE (6 .89) ITER 
WRITE (6.76) 

WRITE(6,79)1P0S 
WRITE  <6( 76 ) 

IMPA1I2  .Ut  1.0)60  TO  7 
WRITE (6.80) 

GO  TO  e 

7  WRITE(6,81) 

8  WRITE (6.76) 

WRITE (6.83) 

WRITE <6.76) 

•0  10  I »1  .ngraph 

10  GRAPHU/wBLANK 

DO  11  1*1 ,NGRAPH,13 

11  6RAPM(I )»DCT 

write (6 ,111 /Graph 

WR1TE(6,111)GRAPH 
DO  13  I«1,NGPAPH 

15  6AAPh<1)«00T 
WRITE <6, 111 /GRAPH 
DO  16  I«1 ,NGRAPH 

16  6RAPh(1)*0LANK 
••••GET  NORMALIZATION 

XRAXwJ.fc 

DO  33  I»NMIN,NMAX 
*AR6*CU(I)*C0NJG(CU(D) 


100 


1 


I  T  (MRS  ,(T.IMX)lKAIMMt 
33  CONTINUE 

XMAX*SGRT(XMAX) 

XCNTR"(MHXN«NHAA>/2 

RHSE  "ATAN (AINAS (CU(ICNTR) ) /REAL (CU(XCNTR))) 

CAAG"CHRLX(w.CfPH$E) 

CIIKX*XMI*CtXR(CAR() 

{••••••SET  «R  MAIN  RIOT  LOOP 

»o  20d  i-nnxn.nmax.incr 

6RARH(N6RARH)"00T 

CAR6"CU(I)/C0ENH 

RHSEMTAN  (AINAS  (CARA) /REAL  (CARS)) 

A0RX"CABS(CAR6> 

IT(RA1I2  .ST.  1.0)AORX>AORI**2 
{•••••••print  an  X 

J"AORI*XNSRMl  O.C 

IMJ  ,0T.  N&RAPH)  J-N6RAPN 

GRARHU)>X 

K>XN6RN1/2.b«RHSE/RX*XNGRH1*1.e 
IT (X  .IT.  1)  K>1 
XT (K  .CT.  N6RARH )  K-NGRARH 
6RARH(K)>RLU$ 

6RARH(1)"00T 
WR1TE(6, 112)1. GRAPH 
GRAPH(J) "BLANK 
6RARH(K)«BLANK 
200  CONTINUE 

00  211  X"1 .NGRRRH 

211  GRARN(l)"0OT 
WRITCCS.IIOGRARM 
00  212  I"1 .NGRARH 

212  GRARN(I)"BLANK 

00  213  I"1.NGRARN,1S 

213  ARAPMO  )«00T 
MRXTCCA.IIOGRARH 
WRITE (4. Ill) GRARN 
RETURN 

CN0 

SUBROUTINE  GAINM2 (U2P.NT0T ,NM1R  ,GLC .QRGLP.6N0FS > 

C ••••SMALL- SIGNAL  GAIN  XS  GlC  AT  CENTER  ANO  REOUCES  OUAORATIC ALLT  TO 
C  GLC«(1.0-006tP)  AT  EOCES (PROVIDED  GNOTS«G.O).  XT  SNOTS  .N(.  0.7. 
C  CENTER  OF  BUA0RAT1C  GAIK  PROFILE  XS  SMXTTEB  BT  THE  TRACTION  SNOTS 
C  OF  THE  MALT-WIDTH  OF  MIRROR  2. 

LEVEL  2.  U2R 
COMPLEX  U2PCNT0T) 

XNMIR-NMIR 

XNTOT-NTOT 

XH1D"(XNT0T«1. 0/2.0 

LXML"(NTCT-NMXR)/2P1 

LIMU"NT0T-LXML*1 

00  29AS  I"LlMLtLlMU 

XT (006LR  .EG.  0.0)60  TO  T 

XXMIO"XNXO-SNOFS*(XNPXR/2.C) 

«x»x 

BR6"((XX-XXHlD>/XNPXR*2)**2 

SL"6LC*(1.0-OR6*OOSLR) 

SO  TO  0 
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7  SL*GLC 

3  BET1*U2P(X)*CONJS(U2P(X>> 

ICNT*1 
BET4*aETl 
1  SAVE"BET4 

BET4*tL«BET1-j.S*AL06<BET4/t>|T1) 

BET4*<BET4*SAVE)/2.£ 

TEST* ABSIBET4 -SAVE) 

TEST  MADE  RELATIVE  ON  1-31-83 
TEST*TEST/SAVE 
XMTEST-C.3uOu1)20,2C,1O 
10  XCNT*1CNT+1 

IMlCNT-50)1T,i;C,10C 

17  IF(BET4)1b9,1O0.1 
100  ICNT»1 

BCT4*BET1 
let  SAV£*t*ET4 

XARG*(6L*BET1-BET4)*2.0 

BET4*BCT1*EXP(XARG> 

BET4*(BET4«SAVE)/«.G 

TEST*ABS(BET4-SAVE) 

te$t*test/$ave 

I  MTEST-O.CUOC/I)  20,20,110 
110  ICNT*ICNT*1 

IHXCnT-S9)1G1( 18,18 

18  VR1TE<6,83)GL.BET1,8ET4 

83  FORHATCIX.ISNPROBlEM  IN  BET4,SriC.S> 

STOP 

20  AIIG*SQRT (BETA) 

C4S  U2P(X)-u2P<I>/CABS<U2P<1))*ARG 
RETURN 
END 

SUBROUTINE  BOUALCNTOT.NHXR.U) 

LEVEL  2,  U 

COMPLEX  U(NTOT) ,CARG 
XNM1N-NM1R 

L1ML*<NT0T-NMIR)/2«1 
L1MU*NT0T-LIML«1 
CARG«CNPLX<0. 0,0.0) 

RARG"C.C 

XINTS-0.9 

DO  1C  I*LXMt.LlMU 

CAR6*CAR6«U(X) 

RARC*RAR6*C SQRT(U(I)*C0NJ6(U(1 )) ) 

XINTS  •  XXNTS  ♦  U ( I ) *C0NJ6 (U ( X  >  > 

10  CONTINUE 

TVft«L*CARG*CONJG(CARG)/RARG**2 

MRXTC(6,*:)TVB«L 

90  PONHAT(1NO,6NTVB«L*»MO.S) 

CAR6*CARC/XNM1R 

XINTS*XINTS/XNM1R 

TVBttT  *  CAR6  •  COhJG(CA*S)/XINTS 
NR ITS (6,91 )TVBDT 

91  PORMATC*  TVBOT*',M0.S) 

RETURN 

END 

SUBROUTINE  MSLDl  (UNPLD«NTOT ,NMXR ,EPSL,XNAC2«A10A2*XBITR, 
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imrr^MMiCrMiiiircSiircirhMCKti) 

C.....TEMPOHAAT  CHANCES  HAVE  KtN  MADE  ON  02-12-01  IN  011(1  TO  (MM 
C*a*aaPR0PA6ATX0N  (AON  ONE  HALT  Of  (ULL  OUTPUT  APERTURE. 

e 

C*aa**NTOTH  IS  NAK  Of  NTOT. 

C 

C«****ANNIAM  IS  MALE  OS  XNRXR. 

C 

. . UPPER  LIMITS  0(  SUMMATION  CHANCED  (RON  NTOT  TO  NTOTh  AT  (R2037  AND 

C***»«rR2d«3>  XNMIRN  REPLACED  XNHIR  AT  (R?w5S«rP2C56«AND  (R2381. 
C*****ALSO  rAKTR  AT  (02022*1  AND  (R2325  AND  (R21C0 
LEVEL  2.  UN(LD 
DIMENSION  UNTLD (NTOT ) 

DIMENSION  IMX(r(2«j0),rrNAX<20G) 

DIMENSION  6R  APM (101  )  a  (flOCI) 

complex  crcNtrrrcT.ArrAH^rrAP'AMCPx 

COMPLEX  UN(lD 
COMPLEX  CrCTRf 
DATA  Pt/3. 141502654/ 

76  PORMAT(lHO) 

77  rORMATdHlX 

111  rORMAT(20X,101A1) 

112  rOAHATMX,  (10.5,  4X,  13,  2X,  10lAl> 

DATA  PI/3. 14150265367 

DATA  BLANK/lM  / , DOT /iH . / , X /1NX / 

N(PLD*101 

XNMXR-NMIR 

XNTOTaNTOT 

(ACTR«XNA62/AlOA2 

0((ST«EPSL*XNMIR/2.0*<(ACTR-1.0>/(ACTA 

XM1D>(VNTOT«1.0)/2.0 

X2«XNID-0((ST-XNMIR/2.0/(ACTR 

X3>RMI0*O(rST«RNNIR/2.0/(ACTR 

NNorr-<N((LO*i)/2 

( AKTRal .0 

IMSPCLrr  .ED.  1.0)rAKTR-2.0/(1.0>A10A2) 

I((SPCL(r  .ED.  2.0) (AKTRB1 .O/ERCMSK 
ircsPCLrr  .eo.  i.cmcTA-i.CMC 
XCNTR  a  XNT0T/2.0 

KtSPCLrr  .16.  1.G>XCMRaXCNTA-XNHIR/4.0«<1.J*1.G/(ACTA> 

XNHX  •  XNHXR 
XNMlRHaXNMI R /£ .0 

ircsPCLrr  .ct.  o.oixami  a  xnmzrm 
limml«i 

LXMMU  a  NTOT 
NTOTM-NTOT/2 

l((SPCLrr  .ED.  l.uKXMPU  a  NTOTH 
IT  (SPCLM  .NE.  2.0)60  TO  777 
XMIDLaXNTOT/2.0  -  XNM1R/4.0*(1.G*A10A£) 

IR6MSK  a  XHNIR*rRCMSK/2.C 
LIMML*  XM1DL  •  IR6MSK 
LIMMUa  XMIDL  ♦  XR6NSK 
XNMI  •  LIMPU  -  LlMNL  ♦  1 
NRITE<6,80>IR6MSK,LIMPL,LIMPU,XAMI 

•9  (0RHAT(3I6»(1t.1> 

777  CONTINUE 

DO  363  1*1  ,NMLD 


r*i  TPaVAts/ 1  o.  0«run 

. . par-field  angle  it  in  units  or  lanbda/vidth.  vidth  is  geometric-* 

c*****optic$  she  or  enitteo  bean,  i,  e.  width  or  last  reflecting******* 

C****«  CLAA6ER  rOA  CONFOCAl  CASE)  PIAROR  TINES  NA6NIFICAT10N 
C****. PROPAGATING  PAST  0BSCURIN6  P1AROA  TO  EOCUS1NC  NX AAOA ( JF  AvY) »••••• 
AR6a-THETP*2.U*PI/XNNlR 
AR6S*S1N(AAG) 

4R6C*C0S  C ARC) 

CFCNaCNPLXtARbC.AAGS) 
fffct-chplx<i. 0,0.0) 
ae  r an  *chplx(c. 0,0.0) 

srrAN  *cnplx(q. c,o.O) 

DO  350  J"LXNNL,LINNU 
XJ«J 

IJaRJ-C.S 

CPCTAE  *  CNPIX (1 .0,0*0) 

IFtFNFCS  .EG.  3.9)60  TO  302 

AA6PH  •  <(XCNTA-XJ>/XNNIR*2.0)**2*PI*FNFCS 

AAGCF  >  COS( A AGPH) 

AR6SF  •  SXNtAAGPM) 

CFCTAF  «  CNPLX(AAGCF.«AAGSF) 

3C2  rrrcT«rrfcT»crcN 
rCTA»0.0 

IF<M>X2*1.  0)323,320,310 

310  IF(XJ-X2)311, 311,312 

311  FCTR*XJ*1.0-X2 
60  TO  318 

312  IF CXJ-X3«1.C)31 3.314,314 

313  FCTA«1.U 
60  TO  318 

314  IF(XJ-X3)31S,320,320 
313  FCTA«x3-XJ 

320  CONTINUE 

318  BFFAN  «6 F FAN  ♦FFFCT*UNFLO<J>*FCTR  •  CFCTAF 
AFFAN  *AF  FAN  «FFFCT*UNFLD tj )*CFCTRF 

330  CONTINUE 

AFFAN-AFFAN/XNNI 

BFFAN-BFFAN/XNNI 

ANCPX«AFFAN-BFFAN 

FFIfl)aREAL<ANCPX)**2«AINA6(ANCPX)**2 
360  CONTINUE 

C* ••••CALCULATE  AVERAGE  NEAR-FIELD  PONEA .•*•*••••••**••**•••**••••**•*•• 

AS«NF«0.0 

BS«NF»0.9 

DO  332  J»LlNNL,LlNNU 

XAR6*REAL(UNFl0(J>)**2«A1HAG(UNFL0<J))**2 

XJ«J 

XJ*XJ-0.5 

1F(XJ-X2«1 .0)330, 330,326 
326  IFtRJ-R 2)321,321, 322 

321  FCTR*XJ*1.0-X2 
60  TO  328 

322  XMA2-X3*1.w>323,324,324 

323  FCTA-I.O 
60  TO  328 

324  IMIJ*X3)323,33C,330 


•% 
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325  rtTKlMi 

329  (SaNMISIkf«Olt«fCT* 

330  CONTINUE 
ASQNF*AXQNF«XARb 

352  CONTINUE 

PAVNF*ASGNF-eS0NF 

pavnf*paxnf/xnmi 

PAVNF*PAVNF/(1.(J-1.3SFACTR) 

DEN-PAVNF* <1.0-1. 0/FACT A )**2 
00  370  1*1,nFFL0 
FFI<I)»FFKI)/OEN 
373  CONTINUE 
ITEA«XHITR 

C*****6ET  NORMALIZATION 
XMAX  •  0*0 
IMXFF(ITER)  ■  1 
00  30  1  *  1  .  NF FLO 

1F(FF1<1)  *GT.  XMAX)  IMXFF(ITER)  »  I 
IF(FFKI)  *GT.  XMAX)  XMAX  •  FFKI) 

30  CONTINUE 

FFMAX(ITER)  *  XMAX 
IFECFFLD  .CT.  2.0)  GO  TO  392 
FF 1NT*C  *C 
DO  386  1*1 tNFFLO 
308  FFINT*FF2NT«FF1(I) 

FF1NT*FFINT/10*0* FAKIR 
FFINT*FFINT«(1.C-1.3/FACTR) 

NGRAPH*101 
XN6AH1*NGAAPH-1 
C**...*. PRINT  HEADING 
WRITE (6*77) 

WRITE (6 .79) ITER 
79  F0RHATC63X,  'ITERATION',  I*) 

WRITE  <6 • 76 ) 

73  FORMAT!  60X,  'FAR-FIELD  INTENSITY') 

WRITE16 t?6) 

DO  10  1*1 «NGRAPH 

10  6RAPH<2)*6LANX 

00  11  1*1 tNGRAPM* 10 

11  GRAPMO  )*D0T 
WRITE (6«111 ) GRAPH 
WRITE(6(111)GRAPM 
DO  15  1 *1 (NGRAPH 

15  GRAPH < 1 ) *D0T 
WRITE(6t111)GRAPH 
00  16  1*1 (NGRAPH 

16  GRAPH! I )*6LANK 
C******SET  UP  PAIN  PLOT  LOOP 

00  2C0  1*1.  nffld 

6RAPH(1)*00T 
GRAPH (NGRA PH )*OOT 
K  •  FFKD/XPAX  •  XN6RP1  ♦  1.0 
IF(K  *LT«  1)  X *1 
!F(K  .GT*  NGRAPH)  K*NGRAPH 
CRAPH(K)  ■  X 

WRITE (6 (112)  771(1),  I,  GRAPH 
GAApH(K.)*BLANK 
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203  CONTINUE 

•C  211  1*1  ,N€AAM 

211  6AAAM(1)*D0T 
HIITI(6,111)6AAAH 
60  212  X*1*N6AAPN 

212  6AAAN(1)*BLANK 

»0  213  l*1,N6AA9M,la 

213  6AAf»H(D*60T 
niTtUillIXlAN 
HAITI (At 111 )6AAPM 
HAITI (6,391)PAVNF 

391  F0AHAT(1X,'PAVNF*',E1S.S) 

HAITI (6*39P)FF1NT 

390  F©ANAT(1X,'FFXNT*'.E1S.S) 

392  CONTINUE 
AC TUAN 
ChO 

5UBAOUTXNC  NFLPP  (I  A*  FF ,  F  MAX  .ANITA  ,  ISAM  , 

1 XNF2  »XNAS1 tEASL* ANNX A t TICCA  f ACNHI > 

01NCNSI0N  ZMXFF(2C0)  •  I  IN  AX  (2  POX 
DIMENSION  6 AAPH( 101 ) 

OATA  HLANK/lH  /.DOT/lH./tX/lHX/,ZEAO/lHw/(Pi.US/1M*/ 

75  FOAMaTC  NOOI-FOANATION  FA A-F I  ELD  AESULTS.') 

76  FOAMATdN  ) 

77  FOAAAT(IHI) 

79  FOAMAT('  XNF2*  '.FtS.S,'  XNA61*  ',110.5,'  EPSl* 

1F10.5,'  XNNIA*  ',F10. 1 • '  TSCED*  ',F2C.0> 

60  FOAMAT(19x,3hO.O,37x,'XKX  IP/XG  XXX',37X,3H1.0> 

61  FOANAT(18x,*H«5.0,37x, '♦♦♦  TNETA  ♦♦♦', 37X.3M5 .0) 

65  FOAHAT(2I9,5F10«5) 

86  FOAHATC'  ITCAATXON  XMXFF  FPNAX(I)') 

111  FOAHAT (2CX,101A1) 

112  FOAMAT(15X«I3«2X,101A1) 

XTCA  *  XHXTA 

HAITI <6,77) 

HAITI <6,7$) 

HAITI <6,76) 

HAlTt(6«66) 

HA1TCC6,76> 

,  60  ICQ  1*1, XTCA 
HAITI (6,85) X,XNXFF(X),FFNAX (I) 

ICO  CONTINUE 
HAITI (6,77) 

HAITI (6, 79)  XNF2 ,XHA6l ,EPSL«XNP.  IA , TIE ED 
HA1TK6.76) 

HA1TI(6,8C) 

NA!TE(6,76> 

NtAAPH*1C1 
XN6AH1*N6AAPM-1 
60  10  I*1,N«AAPH 
1C-6AAPH(t)*BLANK 

60  11  1*1 , NHAAfH ,10 
11  «AAPM(1)*60T 

HAITI(6,111)«AAPN 
HAITI  CO#  111)  66  AM 
60  IS  1*1 fNOAAAN 
is  OAADH(i)  ■  DOT 


WRITE (6,111)GRAPH 
DO  16  1»1 , NGRAPm 
1ft  CRAPH(I)*BLANK 
90  I*1,1T£9 

SRAPH(1)*D0T 
GRAPH(NGRAPH)*90T 
J  •  IHXFF(I) 

IF (J  .GT.  NGRAPH)  J  *  NGRAPH 
GR APN ( J )*PLU$ 

K  •  FFPAX(I)*XNGRH1«1.C 
IKK  .IT  .1)  K  •  1 
IKK  .GT.  NGRAPH)  K  »  NGRAPH 
GRAPH(K)*X 
WRITE (6.112)1 .GRAPH 
GRAPH(J) "BLANK 
GRAPH (K ) *61 ANK 
I  F  ( ISKHF  .EC.  0)  GO  TO  23u 
90  19C  I$«1, ISKHF 
19C  WRITE (6.76) 

23 0  CONTINUE 

00  211  1*1 .NGRAPH 

211  GRAPH(1)*D0T 
WRITE (6,111 )GRAPH 
90  212  1*1, NGRAPH 

212  GRAPH (I }*BLANK 

90  213  I«1, NGRAPH, 10 

213  GRAPH(I)*DOT 
WR1TE(6,111)6RAPH 
WRITE  (6,1 1DGRAPH 
WR ITS (6 ,76 ) 

WR1TE(6,81 > 

1 F (PCMHF  .LE.  0.0)60  TO  21* 

WRITE (1,*1C)XHITA,XNF2,XNAG1,EPSL,XNPIR,YSEED 
*10  F0RPAT(S(F1*.S/),F1*.S) 

90  300  1*1  i  ITER 
RCOUNT  ■  I 
RHXFF  *  IHXFF(l) 

WRITEd ,*00)RC0UNT,RPXFF,FFHAX(I) 

*00  FORHAT (3F10.3 ) 

300  CONTINUE 
21*  RETURN 
SNO 

SUBROUTINE  FOUNT (DATA, NN,N9IN»ISIGN, I  FOR", WORK) 

LEVEL  2,  9ATA 
9ESCR1PTI0N 

FOUNT  PERFORNS  AN  N-DINENSIONAL  FAST  FOURIER  TRANSFORM  ON  AN 
N-DIHENSIONAL  ARRAT  OF  COHPLLX  DATA.  THE  TRANSFQRH  PERF0RNE9 
HAT  BE  EXPRESSED  AS  FOLLCwS  — 

TRANSFORH(J1,02,...)*SUH(OATA(l1,I2,...)*fc1**((I1»1)*(Ol*1)) 

•W***((I2-1)«(J2-D) 

•  .  .  .  ) 

WHERE  H  AN9  01  RUN  FRON  1  TO  NN(1),  AND 
I?  AND  22  RUN  FRO"  1  TO  NN(2),  ETC. 

AND 

fcl  •  EXP(IS!GN*2*PI*SQRT(-1)/Nh(1>)  ,  etc. 
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FOR  ONE  DIMENSION,  The  transform  is  precisely 

TRANSFCRM(jI)  ■  $UH(DATA(11)*w1**((I1-1)*(J1>1))) 

FOuRT  IS  FASTEST  WHEN  THE  NUMBER  OF  DATA  VALUES  IN  EACH 
DIMENSION  IS  A  HIGHLY  COMPOSITE  (FACTORABLE)  NUMBER. 

INPUT 

GIVEN  AN  N-OIMENSIONAL  COMPLEX  ARRAY  OF  DATA 

DATA  —  COMPLEX  ARRAY  IN  WHICH  THE  DATA  TO  BE  TRANSFORMED 

is  placed,  upon  return  to  calling  program  data 

CONTAINS  THE  TRANSFORM  VALUES. 

NN  -  INTEGER  ARRAY  GIVING  THE  (POSITIVE)  NUMBER  OF  POINTS, 
OR  VALUES,  IN  EACH  DIMENSION,  RESPECTIVELY. 

NDIM  -  NUMBER  OF  DIMENSIONS  (INTEGER)  NCIM.GE.1 

1SIGN  -  INTEGER  GIVING  DIRECTION  OF  TRANSFORM  TO  BE  DONE. 

»  -1  IMPLIES  FORWARD 
»  «1  IMPLIES  BACKWARD 

I  FORM  -  INTEGER  PARAMETER  DESCRIBING  THE  FORM  OF  THE  DATA. 

«  1  IMPLIES  THE  DATA  IS  COMPLEX  (NON-TRI VI ALLY) . 

■  C  IMPLIES  THE  DATA  IS  ACTUALLY  REAL.  I.E.,  THE 
IMAGINARY  PART  OF  EACH  COMPLEX  ELEMENT  OF  DATA 
IS  ZERO.  FOURT  IS  SIGNIFICANTLY  FASTER  WHEN 
1F0RH«C. 

WORK  -  COMPLEX  WORK  ARRAY.  WORK  MUST  BE  DIMENSIONED  AS 

LARGE  AS  THE  LARGEST  DIMENSION  OF  DATA  WHICH  IS  NOT 
A  POWER  OF  TWO.  IF  ALL  DIMENSIONS  OF  DATA  ARE 
POWERS  OF  TWO  THEN  WORK  IS  A  DUMMY  ARGUMENT. 

OUTPUT 

TRANSFORMED  VALUES  OF  DATA  ARRAY. IF  COMBINED  FORWARD  AND 
BACKWARD  TRANSFORMS  ARE  MADE  IN  sequence, THE  VALUES  IN 
THE  DATA  ARRAY  ARE  THE  ORIGINAL  VALUES  MULTIPLIED  BY 
NN(1)*NN(2)*  .  .  .  NN(NDIM) 

REMARKS 

control  DATA  6603  VERSION 

CONSULT  REFERENCE  IN  MAKING  CHOICES  OF  TRANSFORMED  VARIABLES 
OR  FREQUENCIES  IN  DOING  CALCULATIONS  IN  THE  TRANSFORMED 
PLANE 
REFERENCE 

THREE  FORTRAN  PROGRAMS  THAT  PERFORM  THE 
COOLEY-TUKEY  FOURIER  TRANSFORM 
N.  M.  BRENNER 

M.  I.  T.  LINCOLN  LABORATORIES 
TECHNICAL  NOTE  19e7-2 
2!  JULY  1967 
AO  657319 

SUBROUTINE  FOURT (D A TA,NN,. NDIM, I  SIGN, I  FORM, GORKI 
DIMENSION  DATA ( 1 ) ,NN (1 >  ,1  FACT (3i  ),  WORK ( 1  ) 

TWOPI  »  6.26J1853U717959 
RTNLF  »  0. 7*710678118655 

1  NTOT«2 

DO  i  IDIM«1,NDIM 

2  NTOT«NW»NN(IDIN) 

NP1«2 

DO  91C  IDIM«1,NDIM 
N«NN(IB1N) 
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*P2»NPl»N 

IF<h. £0.1)60  TO  900 
5  *«N 
NTwO«*P1 
I  F»  1 
IDIV«2 

1?  IQUOT-P/IOIW 

IBEK-M-1DIVM8U0T 
IF(iauOT-IOIV>53,11 ,11 

11  If <IRfcP)2C,12,23 

12  NTWO*NTwO*NTWO 
IFACT<IF)»IDIV 
I F*X  F*1 

"■loyor 
60T0  10 
20  IDIV-3 
IN0N2-IF 

30  IQUOT»P/IDIV 
IREM«H-I01V«I8U0T 
IF(XOUOT-1DXV)60,31  ,31 

31  1F<IREM)40,32,4G 

32  IFACT<IF)»IOIV 
IF«IF*1 
*■10001 

GOTO  30 

40  I0IWIDIV*2 
SOTO  30 

50  I N0N2 ■ X F 
IF<XR£B)60, 51,60 

51  NTUO-NTtaO^NTWO 
GOTO  70 

60  IF»CT<IE)»H 

70  NOh2P*NP2/NTUO 
IC ASE*1 
IFPIN«1 
IlRNG*NPl 

IFtXOXP “4)74, 100,100 
74  IF(IF0RH)71 ,71 ,100 

71  ICASE-2 

I1RNG«NPC»(1«f,PREV/2) 

IF(IDI«-1)72,72,100 

72  IC*$E«3 

I 1 RNG»tPl 

IF(l,TyO-NPl)1w0,1w0,73 

73  IC»SE«4 
IF*XN«2 
NTwO*NThO/2 
N-N/2 
NP2»NP272 
ntot»ntct/2 
I»1 

00  it  4«1 .'•TOT 
DATA(J)»0ATA(1) 

20  I » I ♦ 2 

K3  IF(M0*2P-1)131,101,20o 
1*1  NP2HFaNP2/2 
J»1 
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00  1S«>  12*1,NP2,NP1 
lf(J-I2>l21,13C,13C 
121  I1*AX*12*hP1-2 

DC  125  11*I2,llKAX,2 
00  125  1 3-X 1  ,NT0TtNP2 
J3*J*l3-li 
TENPR  *0ATA  < 13  > 
TE*»PI*0ATA(I3*1> 
0ATA(I3)«0*TA(J3) 
DATA(13«1)*0ATA< J3*1> 
DATA(j3>*T£*PR 
125  0ATACJ3«1)«TEBP1 
133  N*NP2Hf 

U0  tfO-N)1Sw,1Sw,1*1 
1*1  J*J-P 
**»/2 

If («-NP1)150, 1*0.142 
153  J*J*P 

GOTO  3 Cw 
200  NOORK  *2  *N 

00  270  I1*1,NP1,2 
DO  273  13*11 , nTOT |NP2 
J*I3 

DO  26C  I  *  1 » NOORK  1 2 
1FCICASE-3121U. 220,210 
210  *ORK(l)*DATA(J) 

bORK(I«1)*DATA(J«1> 

GOTO  2*C 

220  W0RK(I)*DATA(J) 

UORK(I«1 1*3,0 
2*0  lfP2*NP2 
I f  *1 f  MI M 

250  IfP1-IfP2/IfACT(If) 
J*J«1FP1 

1F<J-I3-1FP2126C,2S5,255 
255  J-J-1FP2 
IfP2«IfP1 
If«If*1 

If (IfP2-NPl)260,26C,25C 
260  CONTINUE 

I2**AX*I2«NP£-NP1 

1*1 

OO  270  i2«I3,i2*ax,np1 

OATA(I2)*bORK(l> 

OATA(12«1)«uORK(I*1) 

273  I»I*2 

300  1MNTwO«NPT)6w0.603,?C5 
3.5  NP1 TU*NP1 *NP1 
IPA««NT*C/NP1 
310  If (IPAR-21353.333.3iC 
32r ' IPAR*1PAR/* 

GOTO  310 

320  00  3«w  11*1.11 RNG.2 

00  3*0  Kl*11.NT0T,NP1Tb 
«2*Xl*NPl 

tippr*oata(x2) 

TEPPI*0ATA(X2«1) 
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DATA(K2)a9ATAU1)~TERPR 
DATAU2«1)aDATAU1*1)-URPX 
DATA  U1>  ■CAT AU1>*TE*P« 

343  DATAO 1«1)«0ATA(K>*1)*TEMPI 
550  lRAX«hM 

35 i  IF <RRAX-hTbC/i)J7w,6:o,tCO 
373  LrAX«MAXC(NPlTUtMAl/L) 

DO  57b  LaNP1,LRAX,NP1Tb 
*aL 

IF (rrax-NPI )*<C,*2C,35C 
3d0  TNETAa-TUOPI ‘FLOAT <L)/ FLOAT U»HR AX) 
IF(ISIGN)43C, 393,393 
393  ThET  Aa“THET  A 
433  WRaCOSCTMETA) 

KIaSIRCTHETA) 

413  b£RaWR‘bR-»II‘RX 
«2Ia2.0‘bR‘Rl 
W3Raw2»*WR-b2I*WX 
u31aU2P*wi*w2X*bR 
42C  DO  33C  1 1 *1  •  1 1 RN6, 2 
KRXSaH‘IPAR‘R 
IF (PRAX-NP 1)430,430,440 
430  KRIN‘11 
440  KDXFaXPAR‘PMAX 
453  KSTEP*4*KD1F 

I F(KSTEP-HTbO)460, 460,530 
463  DO  520  Kl»KRlh,ATOT,KSTEP 
k2>k1«X01F 
K3>K2«KD1F 
K4‘K2«K0IF 

IF (PPAX-NP1 >470,470,460 
47C  UlR>DATA(Kl)«DATA(K2) 

UlIa0ATAU1«1)«DATAU2‘1) 

U2R‘DATA(K3)«DATA(K4) 

U2IaDATA(K3«1)‘DATA(K4«1) 

U3R«DATA(Kl)-D«TA(K2) 

u3XOATA(k1«1)*DATA(r2«1> 

IF (!$XGN>471 ,472,472 

471  U4RaDATAU3‘1)-0ATAU4‘1) 
U4X«0ATA(K4)-DATA(k3) 

GOTO  Slw 

472  U4P*DATA(K4«1 )>DATA(R3«1) 
U4X«DATAU3)-DATA(k4) 

GOTO  51C 

463  T2Raw2R‘DATA<K2)*u2X‘DATAU2«1) 
*2!>V2R*0ATA(k2«1)«u21*DATA(k2) 
T2A‘WA‘DATA(K2)-*1‘6ATAU3*1> 
T3X«4R‘6AT4(k3‘1)*p1«DATAU3> 
T4R>u3R*DATA(K4)«|,3l*6ATA(K4«1) 
T4laW3R*DATA(K4«1)*w!l*DATA(K4) 

UlRaDATA(Kl)‘T2A 

UlIaDATA<*1‘1)‘T2l 

U£Pa?J*‘T4R 

UiIaT3X‘T4I 

U3RaDATA(Xl )»T2R 

U3XaDATA(x1*1)“T2i 

IF(XSXGN>49<,S33,50« 
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490  U4MV3I-T4! 

U41*T*B-T3« 

SOTO  510 
SCO  U4B*T4I«T3l 
II4I*T3*-T4« 

51C  BATAU1>*UlR*U2ft 
6ATaU1*1>*U11«U21 
BaTa(k2>*U3R«U4R 
BATA(K2«1)*u31«U4X 
BATA 

BATaC*3*1>*u1I-U2X 
BATA(K4)*U3R-»U4R 
520  BATA(K4«1)*u31-U41 
MIMKSTIR 
KRIN»4*UniN-Xl>«Xl 
GOTO  4SU 
530  CONTINUE 
N*N*LNA* 

IF<FI-NNAR>549,S40,S70 
540  IMtS16N>550,34C.S63 
550  T£*M*NR 

WR*(MA«I*1>*ATMLF 
Ul*<Wl>TENFft)*RTHLF 
GOTO  41C 
560  UflM-MA 

WP1*(W«-*I>*ATMLF 
NI*<TENPR*nI>*RTHLF 
GOTO  410 
570  CONTINUE 
IPAR*3*IPAR 

GOTO  360 

6UC  1MN0N2P~1  >790,700*601 
601  IFM*NT*0 
IMINON2 

610  1FP2*IFACT<1F>*IFM 

INST  A*»TW0PI /FLOAT (IFACT(IF)) 
IFOSXGN  >612*611*611 

611  TH£TA*-THETA 

612  NSTM*COS<TH£TA> 

NBTPI*SINITH|TA> 

BO  650  41*1  * IFM  »NP1 

THtTN«-TMQPl»FL0AT<4l-1>/Fl0AT<irp2> 

IF ( I SIGN) 614*613 *61 3 

613  TNtTN*-THETN 

614  NNIN«*CO«(TNETN> 

NNINI*SIN(THCTN> 

1  IN  A  11*4 1*1 1RNfc*2 

BO  63C  I1*4l.llNAX,2 
BO  650  13*11, NTOT.NP2 
1*1 

*R*WF>1NR 

MI*NN1NX 

4*NA»*I3*IfP|-irP1 
BO  640  42*13.42NAX,1FM 
TuOhB*MB«hN 
iSNAR^'N^-ir^ 


BO  63:  J3«J2,j3*Aa,irr2 
J«JI»l*MM2-I»M 

SR*PATA(J) 

SI>PATA(J«1) 

OLBS6-0.0 

OlOSI-C.O 

629  StDM'XI 
ST«P1"$J 

SR>TuOwR*SR-OkP$R*P ATA<J) 

SIaTV0MR*S I'OLPS 1*0  ATA ( J ♦ 1 ) 

OLBSR«STRRR 

OCRS I“$TRP| 

J»J-XPP1 

IMJ-JPIM)621,621,62C 
621  wgRK<n>WR*SR-ttl*&I>OLPSft«OATA<J) 

UORK(1«1)*m1*SR«hR*S1-CLPS1«PATA(J*1) 
639  I«X*2 

UT£MMVA»WSTF1 
UR*WR*USTPR-tfX*WS1Pl 
660  WI«*I*WSTPR*«T£RP 
X»1 

00  650  J2*I3,J2*AX,1FP1 
J3*AX»J2*»iP2-lFP2 
00  650  J3»J2, J3NAX.1FP2 
0ATA(j3)*to0RK<l) 

DATA(j3*1)ato0KK(I+1) 

IMIM1 
659  I»I*2 

IFP1»XFP2 

ir(iPpi*NP2i6iof7cotroc 

709  6OTO(90C,aC0,9CC.7C1),ICASE 
7J1  NHALf«N 

N>N«R 

Tnl1Aa-1*0PXm0AHH> 
tf<ISICN)795 (792(792 
7U2  ThETA«-TH£TA 
703  WSTPR-COICTHITA) 

W5TP1«S2MTN|TA) 

MR*mSTPR 

«*I«rSTPI 

IRINaS 

J*IN«2«*HAl.F-1 
60TC  723 

710  J»J*X* 

60  720  1«IPXR,«iTOT,NP* 

3U6R*(BATA(1)*0ATA( J))72*0 

tUPI*10ATA( 1*1 )*0A1A1 J*1 ) ) /2 ( 0 

017B"l0A7AtI >*0AtAl J l >72 • J 

6j»l«(PAlAU*1>-0ATAO*1>)/2.a 

T|APMPI*>PM««l*tlM 

T(MI<kt*iUM*kR*IIH 

ttlAUI'lkPRMIAPl 

0ATA<J«1>*01Fl*T|ftPX 

0I1«U)*I«PI>TIWI 

OAT«7»IImO!M«T|»PI 
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723  J*1*np2 

IMIN*1PX*«2 
J»W»*JA1»i-2 
tfT(RP*W**4$TPI 
WR*wll*tfSTP*-Wl*wSTPl 
Vl*wX*MSTM*UTENP 
72 5  IMIM1N*JN1N>71C.73C.740 
733  IM1$16N>731, 743,740 
731  00  73S  1*IR1N,NT0T  «hPi 
73$  0ATA<X4l>**0ATA(14l> 

743  «P2-NP24NP2 

4T0T*KT0T4*T0T 
J*NT0T4l 
XMAX«NT0T/2«1 
74$  IMN*XRAX-24NNALf 
X-1NXM 
60T0  7SS 

750  0ATA(J>*0ATA(1> 

OATA(J«1)a-OAtA<X«1> 

7S$  I»X42 
J*J-2 

IM1-IMAX)7SC, 760,760 

760  0ATA(J>*0ATAflRI*>-0ATA(lRI»i4l> 
*  »ATA(J«1)-3.C 
XMI-J)770, 780,780 
76$  DATACJ)-DATA(I) 

DATA(J«1)-»ATA(1«1) 

77P  1*1-2 
l*J-2 

If  <1-1*11077$, 775, 7*5 
77$  DATAU)>0ATAflMlNI«BATAUNlN«1> 
DATAU*1)«C.O 

xnAxai*m 

60T0  74$ 

730  DATA(1)*DATA(1)«BATA<2> 
DATA12I-0.3 
GOTO  9CC 

300  IMXl«AG-NPl>805,900,900 
80$  00  340  13*1,NT0T,MP2 
I2AAXaI3*XP2-NP1 
00  863  12*13. I2HAX.NP1 
IMI**I2«I1*NC 
IIUX*12«NM*2 
JHAX*2*XS«NP1-XMXN 
Xf (12*13 >82 0.62 0.8 1C 

810  JRAX*JHAX4ftp2 

82C  If (I 0IP*2 >850,850 ,830 

830  J*3RAX«NPw 

00  840  X*I*1*,XMAX,2 
OAT A ( 1 >*0ATA(1> 
'9ATA(X4l>*-0*TAO«1> 

34?  1*3-2 
333  1*1RAX 

00  84C  X*1PI A.XAAX.APw 

0ATA(X)*0ATA(1) 

0ATaI1*1>*-»ATAN«1> 

863  2*1*R0C 


LIST  0?  SYMBOLS 


8  f  8X »  8j 

81  >  82 


X*  V  X’’  V 


active  Midi*  dimension 
feedback  mirror  half  dimension 
■agnatic  induction  vector 
epeed  of  ll£it  in  vacuus 
electric  displacement  vector 
laser  aperture  output  dimension 
electric  field  vector 
resonator  waist  scale  parameter 
scaled  resonator  focus  distance 
fluorescence 

equivalent  Fresnel  number 
parabola  focal  length 
feedback  mirror  focal  length 
primary  mirror  focal  length 
resonator  "g"  parameters 
statistical  weights 
normalised  absorption  llneshape 
magnetic  field  vector 
Hermlte  polynomial  of  order  ■ 
Hermits  polynomial  of  order  1 
intensity 
Huygens  kernel 
propagation  number 
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LIST  OF  SYMBOLS  (Continued) 


resonator  airror  separation 
active  media  gain  length 
resonator  magnification 
mesh  point  sire 

feedback  mirror  (unless  otherwise  stated) 
primary  mirror  (unless  otherwise  stated) 
tube  Fresnel  number 
lower  laser  level  population 
upper  laser  level  population 
Anan'ev  round  trip  number 
diffraction  limited  round  trip  number 
Fourier  spatial  frequency  component  in  x  coordinate 
Fourier  spatial  frequency  component  in  y  coordinate 
feedback  mirror  radius  of  curvature 
primary  mirror  radius  of  curvature 
resonator  geometric  outcoupllng 
Fourier  transform  in  x  coordinate 
inverse  Fourier  transform  in  x  and  y  coordinates 
spontaneous  emission  lifetime 
lowest  loss  eigenfunction 
U  ^  resonator  eigenfunction 

beam  waist  of  symeietric  confocal  resonator 
beam  waist  of  a  stable  resonator 
laser  beam  dimension 
small  signal  gain 
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LIST  OF  SYMBOLS  (Continued) 


saturable  gain 
Kronecker  delta 
electric  permittivity 
index  of  refraction 

round  trip  angular  radiation  distribution 

lowest  loss  eigenvalue 

complex  eigenvalue 

wavelength  of  light 

magnetic  permeability 

photon  frequency 

time  after  threshold 

complex  amplitude 

transformed  complex  amplitude 

radial  frequency 
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